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We present a calculation of the Nucleon and Delta excited state spectrum on dynamical anisotropic 
clover lattices. A method for operator construction is introduced that allows for the reliable iden- 
tification of the continuum spins of baryon states, overcoming the reduced symmetry of the cubic 
lattice. Using this method, we are able to determine a spectrum of single-particle states for spins 
up to and including J = |, of both parities, the first time this has been achieved in a lattice cal- 
culation. We find a spectrum of states identifiable as admixtures of 5(7(6) ® 0(3) representations 
and a counting of levels that is consistent with the non-relativistic qqq constituent quark model. 
This dense spectrum is incompatible with quark-diquark model solutions to the "missing resonance 
problem" and shows no signs of parity doubling of states. 



I. INTRODUCTION 

Explaining the excitation spectrum of baryons is core 
to our understanding of QCD in the low-energy regime, 
and if we truly understand QCD in the strong-coupling 
regime, we should be able to confront experimental spec- 
troscopic data with first-principles calculations within 
QCD. The experimental investigation of the excited 
baryon spectrum has been a long-standing element of 
the hadronic-physics program. An important goal has 
been the search for so-called "missing resonances" , bary- 
onic states predicted by the quark model based on three 
constituent quarks but which have not yet been ob- 
served experimentally; should such states not be found, 
it may indicate that the baryon spectrum can be mod- 
eled with fewer effective degrees of freedom, such as in 
quark-diquark models. In the past decade, there has 
been an extensive program to collect data on electro- 
magnetic production of one and two mesons at Jefferson 
Lab, MIT-Bates, LEGS, MAMI, ELSA, and GRAAL. To 
analyse these data, and thereby refine our knowledge of 
the baryon spectrum, a variety of physics analysis models 
have been developed at Bonn, George Washington Uni- 
versity, Jefferson Laboratory and Mainz. 

The experimental efforts outlined above should be 
complemented by high-quality ab initio computations 
within lattice QCD. Historically, the calculation of the 
masses of the lowest-lying states, for both baryons and 
mesons, has been a benchmark calculation of this dis- 
cretized, finite-volume computational approach, where 
the aim is well-understood control over the various sys- 
tematic errors that enter into a calculation; for a recent 
review, see [1]. However, there is now increasing effort 



'Electronic address: edwards@jlab.org 

^Electronic address: dudek@jlab.org 

tElectronic address: dgr@jlab.org 

§ Electronic address: stevewal@umd.edu 



aimed at calculating the excited states of the theory, with 
several groups presenting investigations of the low-lying 
excited baryon spectrum, using a variety of discretiza- 
tions, numbers of quark flavors, interpolating operators, 
and fitting methodologies [2HS]. Some aspects of these 
calculations remain unresolved and are the subject of in- 
tense effort, notably the ordering of the Roper resonance 
in the low-lying Nucleon spectrum. 

A basis of baryon operators for states at rest, re- 
specting the (cubic) symmetry of the lattice, was de- 
veloped in Refs. [5J [7], and subsequently used in cal- 
culations of the excited state Nucleon spectrum in 
both quenched QCD 8J, and with two dynamical quark 
flavors [9 . In parallel, we studied Clover fermions on 
anisotropic lattices jTD], [IT], with a finer temporal than 
spatial resolution, enabling the hadron correlation func- 
tions to be observed at short temporal distances and 
hence many energy levels to be extracted. Crucial to 
our determination of the spectrum has been the use of 
the variational method [TMTi] with a large number of in- 
terpolating operators at both the source and the sink; we 
developed and used the "distillation" method, enabling 
the necessary correlation functions to be computed in an 
efficient manner. A recent calculation of the Nucleon, A 
and Q excited-state spectrum demonstrated the efficacy 
of the method[T5]. 

In this paper, we expand the above program of 
computations considerably, extending to baryons the 
spin-identification techniques developed for mesons in 
Refs. [HH [17]. We develop a new basis of interpolat- 
ing operators with good total angular momentum, J, in 
the continuum, which are then subduced to the various 
lattice irreducible representations (irreps). We find that 
the subduced operators retain a memory of their contin- 
uum antecedents to a remarkable degree. For example, 
hadron correlation functions between operators subduced 
from different continuum spins J are suppressed relative 
to those subduced from the same J, illustrating an ap- 
proximate realization of rotational symmetry at the scale 
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TABLE I: The lattice data sets and propagators used in this 
paper. The light and strange quark mass as well as the Q 
baryon mass in temporal lattice units are shown. The pion 
mass, lattice size and number of configurations are listed, as 
well as the number of time-sources and the number of distil- 
lation vectors N VBCS . 



of hadrons. This allows us to determine reliably the spins 
of most of the (single-particle) states, which helps to de- 
lineate between the nearly degenerate energy levels we 
observe in the spectrum. We are thereby able to deter- 
mine the highly excited spectrum, including spins up to 
J = |, and to resolve the masses of the low-lying states 
with a statistical precision at or below 1%. 

The remainder of the paper is organized as follows. In 
Section|n]we describe the parameters of the lattice gauge 
fields used in our calculation. The "distillation" method, 
and its application to the construction of baryon cor- 
relation functions is outlined in Section IIIII Section IIVI 
describes the procedure for constructing baryon inter- 
polating operators with good continuum spin. Angular 
dependence transforming like orbital angular momenta 
is incorporated through covariant derivatives that trans- 
form as representations with L = 1 and L = 2; a de- 
tailed construction is provided in Appendix |A") We also 
develop subduction matrices that allow the continuum 
operators to be subduced to irreducible representations 
of the cubic group; a derivation for half-integer spins up 
to | is given in Appendix \B\ where tables of the sub- 



duction matrices are also given. Our implementation of 
the variational method is presented in Section |VJ and 
our procedure for determining the spins of the result- 
ing lattice states is described in Section |VI[ which also 
shows tests of the approximate rotational invariance in 
the spectra. The stability of the resultant spectra with 
respect to changes in the analysis method is discussed in 



Section VII We present our results in Section VIII be- 



ginning with a detailed discussion of the spectrum at the 
heaviest of our quark masses before examining the quark 
mass dependence. Discussion of multi-particle states is 



given in Section IX We summarize our findings and pro- 
vide our onclusions in Section Ixl 



state. Hence the signal-to-noise ratio for excited-state 
correlators exhibits increasingly rapid degradation with 
Euclidean time with increasing energy. To ameliorate 
this problem we have adopted a dynamical anisotropic 
lattice formulation whereby the temporal extent is dis- 
cretized with a finer lattice spacing than in the spatial 
directions; this approach avoids the computational cost 
that would come from reducing the spacing in all direc- 
tions, and is core to our excited-state spectroscopy pro- 
gram. Improved gauge and fermion actions are used, 
with two mass-degenerate light dynamical quarks and 
one strange dynamical quark, of masses mi and m s re- 
spectively. Details of the formulation of the actions as 
well as the techniques used to determine the anisotropy 
parameters can be found in Refs. [TUI [TTj . 

The lattices have a spatial lattice spacing a s ~ 
0.123 fm with a temporal lattice spacing approximately 
3.5 times smaller, corresponding to a temporal scale 
a^ 1 ~ 5.6 GeV. In this work, results are presented for 
the light quark baryon spectrum at quark mass param- 
eters a t mi = (-0.0808, -0.0830, -0.0840) and a t m s = 
—0.0743, and lattice sizes of 16 3 x 128 with spatial extent 
r~> 2fm. The bare strange quark mass is held fixed to its 
tuned value of a t m s = —0.0743; some details of the lat- 
tices are provided in Table [TJ The lattice scale, as quoted 
above, is determined by an extrapolation to the physical 
quark mass limit using the Q baryon mass (denoted by 
Q't'm-n)- To facilitate comparisons of the spectrum at dif- 
ferent quark masses, the ratio of hadron masses with the 
n baryon mass obtained on the same ensemble is used to 
remove the explicit scale dependence, following Ref. [llj . 



III. CORRELATOR CONSTRUCTION USING 
DISTILLATION 

The determination of the excited baryon spectrum pro- 
ceeds through the calculation of matrices of correlation 
functions between baryon creation and annihilation op- 
erators at time and t respectively: 

C ij (t) = (0\O i (t)O}(0)\0). 

Inserting a complete set of eigenstates of the Hamilto- 
nian, we have 



C ij (t)=J2^( \ O i\ n )(*\ O ]\ )' 



-E„t 



(1) 



II. GAUGE FIELDS 

A major challenge in the extraction of the spectrum 
of excited states from a lattice calculation is that the 
correlation functions, or more specifically the principal 
correlators of the variational method that correspond 
to excited states, decay increasingly rapidly with Eu- 
clidean time as the energy of the state increases, whilst 
the noise behaves in the same manner as in the ground 



where the sum is over all states that have the same quan- 
tum numbers as the interpolating operators {Oi}. Note 
that in a finite volume, this yields a discrete set of ener- 
gies. 

Smearing is a well-established means of suppressing 
the short-distance modes that do not contribute signif- 
icantly to the low-energy part of the spectrum, and in 
turn, allows for the construction of operators that cou- 
ple predominantly to the low-lying states. A widely 
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adopted version is Jacobi smearing, which uses the three- 
dimensional Laplacian, 



where 
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3 = 1 



where the gauge fields, U may be constructed from an ap- 
propriate covariant gauge-field-smearing algorithm |18j . 
From this a simple smearing operator, 



■W*)= i 



is subsequently applied to the quark fields ip; for 
large n CT , this approaches a Gaussian form exp<rV 2 (i). 
"Distillation" [19] , the method we adopt, replaces the 
smearing function by an outer product over the low-lying 
eigenmodes of the Laplacian, 



N 



□xv(*) = E^* ) (*)e^ t (<), 



(2) 



k=l 



where the ^ (*) is the k th eigenvector of ViL(i), or- 
dered by the magnitude of the eigenvalue; the (volume- 
dependent) number of modes N should be sufficient to 
sample the required angular structure at the hadronic 
scale jTTJ [T5], but is small compared to the number of 
sites on a time slice. Thus distillation is a highly efficient 
way of computing hadron correlation functions. 

To illustrate how distillation is applied to the construc- 
tion of the baryon correlators, we specialize to the case 
of a positively charged isospin-^ baryon. A generic anni- 
hilation operator can be written 

Oi(t) = ^St i ^(T 1 nd)%(T 2 nu) b fi 0? 3 nu) c J (t), (3) 

where u and d are u— and d— quark fields respectively, 
Tj is a spatial operator, including possible displacements, 
acting on quark j, a, b, c are color indices, and a, j3, 7 are 
spin indices; S encodes the spin structure of the operator, 
and is constructed so that the operator has the desired 
quantum numbers, as discussed in the next section. We 
now construct a baryon correlation function as 



Cy(t) = 



abcabc ryi 



q*J_ 
'a-Pl a/37 



[{T^dT^Uuf^Uuf^t)} 

•[(Jnr 1 )^(snr 2 )|(unr 3 )^(o) 



where the bar over S and T indicate these belong to the 
creation operator. Inserting the outer-product decompo- 
sition of □ from Eq. [2] we can express the correlation 
function as 



Ca(t) 



<.<P*r> (t) ^WW (0) 



a/3 7 



^(*,0)r«|(i,0)r^(*,0) 



tffoKT) = e0 6c 5 ^ (riC ( P ) )a(r2eW)6(r ^ 
encodes the choice of operator and 



(r)\c 



is the operator-independent "perambulator" , with p, p 
the eigenvector indices, and M the usual discretized 
Dirac operator. The perambulators play the role of the 
quark propagators between smeared sources and sinks. 
Once the perambulators have been computed, the fac- 
torization exhibited in Eq. Q enables the correlators to 
be computed between any operators expressed through 
including those with displaced quark fields, at 
both the source and the sink. This feature will be key 
to our use of the variational method and the subsequent 
extraction of a spin-identified baryon spectrum. 



IV. CONSTRUCTION OF BARYON 
OPERATORS 



The construction of a comprehensive basis of baryon 
interpolation operators is critical to the successful appli- 
cation of the variational method. The lowest-lying states 
in the spectrum can be captured with color-singlet local 
baryon operators, of the form given in Eq. S 1 . Combi- 
nations of such fields can be formed so as to have definite 
quantum numbers, including definite symmetry proper- 
ties on the cubic lattice. The introduction of a rotation- 
ally symmetric smearing operation, be it Jacobi smearing 
or distillation, in which the quark fields ip(x) are replaced 
by quasi-local, smeared fields □^(x), does not change 
the symmetry properties of the interpolating operators. 
However, such local and quasi-local operators can only 
provide access to states with spins up to |. In order to 
access higher spins in the spectrum, and to have a suffi- 
cient basis of operators to effectively apply the variational 
method, it is necessary to employ non-local operators. 

A basis of baryon interpolating operators incorporat- 
ing angular structure, respecting the cubic symmetry of 
the lattice, and able to access higher spins in the spec- 
trum, was constructed in refs. [HIE], and these operators 
were employed in our earlier determinations of the spec- 
trum. The identification of the continuum spins corre- 
sponding to the calculated energy levels remained chal- 
lenging, however. To overcome this challenge, we adopt 
a different procedure for operator construction: we first 
derive a basis of operators in the continuum, with well- 
defined continuum spin quantum numbers, and then form 
the subduction of those operators into the irreducible rep- 
resentations of the octahedral group of the lattice. 



5 0)rj^(t, 0)r^ (f,0) , (4) 1 

With, in this case, T being just ordinary Dirac gamma matrices. 
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A. Continuum baryon interpolating operators 



all symmetric projector: 



Baryons are color-singlet objects, and thus they in- 
volve totally anti-symmetric combinations of the color 
indices of the three valence quarks. Furthermore baryon 
interpolating operators have to be anti-symmetric un- 
der the exchange of any pair of quarks, which is auto- 
matically satisfied since they are constructed from anti- 
commuting Grassmann fields. Thus the remaining quark 
labels, namely those of flavor, spin and spatial structure, 
have to be in totally symmetric combinations. 

We will construct our baryon interpolating operators 
from products of three quark fields. Before proceeding to 
classify our operators, we note that three objects {cc, y, z} 
can exist in four definite symmetry combinations: sym- 
metric (S), mixed-symmetric (MS), mixed antisymmetric 
(MA) and totally antisymmetric (A) 2 ; projection oper- 
ators that generate these combinations are specified in 
the Appendix, in Eqs. (Al) and (|A2|). We will write 



our baryon interpolating operators by applying projec- 
tion operators that act on the flavor, spin and spatial 
labels of a generic three-quark operator, ipxip2ip3,'- 

B = (Jk F ®«Sr s ^Th^i^iiMz}, 

where J 7 , S and T> are flavor, Dirac spin and spatial pro- 
jection operators, respectively, and the subscripts Zp, Zs 
and specify the symmetry combinations of flavor, 
Dirac spin and spatial labels. For each operator B, we 
must combine the symmetry projection operators such 
that the resulting baryon operator is overall symmetric. 
The rules for combining symmetries of such direct prod- 



ucts are given in Eq. ( A3 1 



To illustrate the construction, we specialize to the case 
of local or quasi-local operators, where the spatial depen- 
dence is the same for each quark, and thereby symmetric. 
Thus we write this simplified interpolating operator as 

B= (Jtf®^)^^}- 

Furthermore, we will only consider the case of two- 
component Pauli spin rather than four-component Dirac 
spin. In the standard convention, Pauli spin involves a 
label s that can take two values: + and — . It is straight- 
forward to extend the construction to Dirac spins using 
the Dirac-Pauli representation of Dirac matrices, as in 
Refs. H>J|7)- Dirac spins involve direct products of two 
Pauli spin representations: one for s-spin and the other 
for p-spin (intrinsic parity). The operator construction 
for Dirac spin is described in Appendix |A") 

The product rules of Eq. (A3) specify three ways to 



combine these flavor and spin projectors to yield an over- 



FaSa 

^75 (-^MS^MS + -^ma^ma) 



(5) 
(6) 
(7) 



The four symmetric spin combinations 5s simply cor- 
respond to the four states or operators of spin |: 

{+ + +}s- {+ + -}s, {+ - -}s, { }s, while the two 

mixed symmetric and two mixed antisymmetric combi- 
nations each correspond to states or operators of spin i; 
there is no antisymmetric combination of three objects 
taking only two values. 

For the case of SU(3) flavor, the flavor-symmetric com- 
bination T$ yields the decuplet (10). Hence Eq. (pSl), with 
J-"s {uuu}, specifies the operators for the spin-^ A ++ . 
The mixed-symmetric combinations J-ms.ma specify the 
octet (8), so that {u(1u}ma,ms correspond to the Nucleon. 
Thus we see that Eq. ^ specifies the operators for the 
spin-| octet. Since there is no 5a spin combination, this 
example cannot provide flavor-singlet (1) interpolating 
operators - angular structure through non-local behavior 
is required. 

Covariant derivatives defined as in Ref. [6] are incorpo- 
rated into the three-quark operators in order to obtain 
suitable representations that transform like orbital an- 
gular momentum. This is necessary in order to obtain 
operators with total angular momentum J > |. First 
the covariant derivatives are combined in definite sym- 
metries with respect to their action on the three quark 
fields. A single derivative is constructed from {D 1 1}e d . 
There are two relevant symmetry combinations, 



L = 1 



D [i] 

^MA 



7e( 2D(3 



D (2)) 



(8) 



where the notation means that the derivative acts 
on the q-th quark. There is no totally antisymmetric 
construction of one derivative and the symmetric combi- 
nation + + is a total derivative that gives 
zero when applied to a baryon with zero momentum - it 
is omitted. Each derivative additionally has a direction 
index (suppressed above) for which we choose a circu- 
lar basis such that they transform under rotations like 
components of a spin-1 object: 

An=±l = ±\(D x ±lD y ) 

D m= o = ~-^D z . 

Totally symmetric baryon operators with one deriva- 
tive are constructed by applying Eq. ^ to the mixed 
symmetry spin-flavor operators according to the rules 
given in Eq. (A3 1, 



2 Mixed symmetry combinations are of definite symmetry under 
exchange of the first two objects. 



y J S ^ I MS^MS ' MAtMA/' 

where superscripts in brackets indicate the number of 
derivatives in the operator. 
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J irreps, A(dim) 



Gi(2) 

H(4)0G 2 (2) 
Gi(2)©ff(4)©G 2 (2) 

d(2)e 1 -ff(4)e 2 #(4) 



TABLE II: Continuum spins subduced into lattice irreps 
A(dim). There are two embeddings of H in J = ~. 



Two derivative operators can be formed in definite 
three-quark symmetry combinations that transform like 
L = 0, 1, 2 as follows, 



L = 


0,2 








, £,[1] £,[1] ) 
MA MA/' 


(9) 


L = 


0,2 






_£,[!] £,[1] 


■ £,[!] £,[1] ) 

MA MA/' 


(10) 


L = 


0,2 






f D [1] D [1] 

' MS MA 


MA MS/' 


(11) 


L 


= 1 






i nW nPJ 

^^MS^MA 


nW nW \ 

^MA^MS^' 


(12) 



The projection into L — 0,1,2 comes from combining 
a pair of derivatives via their (suppressed) direction in- 
dices using an 50(3) Clebsch-Gordan coefficient, i.e. as 
(l,mi; l,m 2 \L, M)D mi D m2 . 

[21 

Although it is allowed, we omit the D s combination 
that couples to L = 0. It corresponds to the spatial 
Laplacian and vanishes at zero momentum. Several pos- 
sibilities occur when totally symmetric baryon operators 
are formed using the rules of Eq. ( A3 ) to combine the 
spatial derivatives and spins, 



4 2] 



[o] 



1 

71' 



(r> [2] ib l0] + n m V; lU| ^ 

l^MSWlS + ^MA^MAj 



[0] 



where no total derivatives are formed from these con- 
structions. The angular momenta of spinors and deriva- 
tives are combined using the standard Clebsch-Gordan 
formula of SU (2) in order to obtain operators with good 
J in the continuum, 

|[J,M])= ^2 | [Ji,mi])(g)| [J 2 ,m 2 ])( JiTOi; J2m 2 \ JM). 
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TABLE III: (left) Number of distinct operators, categorized 
by SU(3)f in each lattice irrep A g>u , using operators with 
up to two derivatives, (right) Number of distinct operators, 
categorized by SU(2)f isospin (I) and strangeness (S) in each 
lattice irrep A 9jU , using operators with up to two derivatives. 
Each operator has dim(A) "rows". There are equal numbers 
of operators in positive and negative parity. 



where subscripts show that the flavor construction is 
mixed symmetry (M), the spin construction also is mixed 
symmetry and the two derivatives are in a symmetric 
combination (S) as in Eq. Direct products of these 
flavor, spin and space constructions yield an overall sym- 
metric set of flavor, spin and spatial labels, as required. 

A spin state like f| ^ can be constructed in several 

ways because Dirac spinors are used and they have four 
components, two upper and two lower components. A 
superscript 1 attached to the spin part indicates that 
the operator is the first of several embeddings with the 
same quantum numbers. Because of this, the basis set 
presented here is over-complete since operators featur- 
ing both derivatives and lower components are specified. 
However, this redundancy is intentional as we need a 
basis with multiple operators in each irrep for use with 
the variational method. If only Pauli spinors were used, 
we would have an SU(6) ® 0(3) classification of opera- 
tors. That subset of our operators is referred to as "non- 
relativistic" . 

We reiterate that the operator construction described 
here provides a basis set for spectrum calculations. We 
are not forcing the symmetries manifested in the oper- 
ator basis into the spectrum results to be shown later. 
Rather, the dynamics of QCD will decide on the eigen- 
states, which can correspond to any linear superpositions 
of the operator basis which may exhibit quite different 
symmetries to those diagonalised by this basis. 



The scheme outlined in this section provides a classi- 
fication of all baryon operators constructed from three 
quarks with either light or strange masses, up to two 
derivatives in their spatial structure and respecting clas- 
sical continuum symmetries. 

A naming convention for such operators is useful. A 
Nucleon operator, with spin and parity of the three 
quarks equal to | , two derivatives coupled into L = 2 



and total spin and parity J 



is denoted as 



(^Vm ® ( 



'm 



D 



[2] 



L=2,S 



(13) 



B. Subduction 

In lattice QCD calculations the theory is discretized 
on a four-dimensional Euclidean grid. The full three- 
dimensional rotational symmetry that classifies energy 
eigenstates in the continuum is reduced to the symmetry 
group of a cube (the cubic symmetry group, or equiv- 
alently the octahedral group). Instead of the infinite 
number of irreducible representations labelled by spin 
and parity, J p , the double-cover cubic group for half- 
integer spin has only six irreducible representations (ir- 



G 



reps): G\ g , H g , G 2g , G lu , H u , G 2u , where parity is de- 
noted by the g (gerade) subscript for positive parity or 
the u (ungerade) subscript for negative parity. The dis- 
tribution of the various M components of a spin- J baryon 
into the lattice irreps is known as subduction, the result 
of which is shown in Table HU 

Extending the analysis of Refs. [THl HZ] to baryons, 
we use "subduction" matrices for half-integer spins to 
project the continuum-based operators to irreducible rep- 
resentations of the octahedral group, 



oli r = J2 s 'C° lJ ' M] > 



(14) 



M 



where 0[ J ' M 1 is a lattice operator constructed as outlined 
in Appendix (A3). As seen in Table [Tij subduction of 
J = I yields two occurrences of the H irrep. Multiple 
occurrences of lattice irreps are the general rule for spins 
higher than 4. We denote by "A the n th occurrence of 
irrep A in the subduction of spin J. For example, the 
two occurrences of octahedral irrep H in the subduction 
of J = | are denoted as and 2 H . For each J — >"A 
there is a subduction matrix S in the values of M and 
the rows of the irrep, r, that maps the continuum spin 
operators to irreducible representations of the octahedral 
group. These subduction matrices are derived for half- 
integer spins up to | in Appendix [Bl and the resulting 
matrices are tabulated in Appendix \a\ 

Classifying operators by SU(3)p symmetry, we have 
the total number of operators in each lattice irrep as 
shown in Table |Hl| lcft). Considered as a broken sym- 
metry with an unbroken SU (2) p isospin symmetry, we 
have the number of operators shown in Table IlFright). 



The operator basis used in this work is constructed using 
SU(3)p flavor symmetry. Of course, the SU(3)p symme- 
try is broken in QCD, most notably by the mass of the 
strange quark. The SU(3) p operator basis allows the fla- 
vor symmetry breaking to be determined by the relative 
degree of overlap onto operators transforming in different 
representations of SU(3)p. This is comparable to what is 
done in quark models and other approaches, which helps 
to relate the lattice results to phenomenology. 

With the construction of operators and their subduc- 
tion to irreps of the cubic group outlined as above, we 
defer further details to the appendices. Appendix [X] pro- 
vides details of the operator constructions, Appendix [B] 
derives the subduction matrices through an analysis of 
the quantum mechanics of continuum spin in the octa- 
hedral representation, and Appendix [B| provides the sub- 
duction matrices for half-integer spins. 



V. CORRELATOR ANALYSIS 

The variational method for spectral extraction [T2T - 
[14], which takes advantage of the multiplicity of op- 
erators within a given symmetry channel to find the 
best (in a variational sense) linear combination of op- 
erators for each state in the spectrum, is now in common 
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FIG. 1: Principal correlator fits according to Eq. (17 1. For 



ten states of the Nucleon H u irrep for ro, = 524 MeV, the 
plots show \"(t) ■ e m " (*-*!)) data and the fit for to = 9. Data 
used in the fit are shown in black, while points excluded from 
the fit are in grey. 



usage[15, 20-23]. Our application of the method follows 
that developed in Refs. [HI [T71 [53] , and applied to the 
analysis of the excited meson spectrum, and we summa- 
rize it here. 

The starting point is the system of generalized eigen- 
value equations for the correlation matrix: 



C(t)v n (t) = X n (t)C(t )v n (t) 



(15) 



where A n (to) = 1, and where there is an orthogonality 
condition on the eigenvectors of different states (n, n'), 
v n — f^n.n' • This orthogonality condition pro- 
vides eigenvectors that distinguish clearly between nearly 
degenerate states, which would be difficult to distinguish 
by their time dependence alone. Equation ( 15 ) is solved 
for eigenvalues, A n , and eigenvectors, v n , independently 
on each timeslice, t. Rather than ordering the states by 
the size of their eigenvalue, which can become uncertain 
because of the high level of degeneracy in the baryon 
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spectrum, we associate states at different timeslices us- 
ing the similarity of their eigenvectors. We choose a ref- 
erence timeslice on which reference eigenvectors are de- 
fined, u" ef = u tt (t re f), and compare eigenvectors on other 

timeslices by finding the maximum value of v*J C(to)v n 
which associates a state n with a reference state n'. Us- 
ing this procedure we observe essentially no "flipping" 
between states in either the principal correlators, A n (i), 
or the eigenvectors, v n (t), as functions of t. 

Any two-point correlation function on a finite spatial 
lattice can be expressed as a spectral decomposition 



C l3 (t) 



y- A Z j e - TOn t 
^ 2m n 



(16) 



where we assume that ( < L t , the temporal length 
of the box, so that the opposite-parity contributions 
arising from the other time ordering on the periodic 
lattice can be ignored. The "overlap factors", Z" = 
(n|C| |0) are related to the eigenvectors by Zf — 

We obtain the masses from fitting the principal cor- 
relators, X n (t), which for large times should tend to 
e -m n (i-t )_ j n p rac tice we allow a second exponential 
in the fit form, and our fit function is 

\ n (t) = (1 - A^e-" 1 '^- 1 ^ + A»e- m '«( t - i »), (17) 

where the fit parameters are ra n , m' n and A n . Typical fits 
for a set of excited states within an irrep are shown in 
Figure [l] where we plot the principal correlator with the 
dominant time-dependence due to state n divided out. If 
a single exponential were to dominate the fit, such a plot 
would show a constant value of unity for all times. For 



the form of Eq. 17 the data would approach a constant 
1 — A at large times, and this is clearly satisfied for t > to- 
Empirically we find that the size of the second expo- 
nential term decreases rapidly as one increases to- Fur- 
ther we find, in agreement with the perturbative analysis 
of Ref. [13] and with our earlier meson analysis, that 
for large to values the m' n extracted are larger than the 
value of TO n=< jim(C) ) the largest "first" exponential mass 
extracted. At smaller t values this is not necessarily 
true and is indicative of an insufficient number of states 



in Eq. (161 to describe C(t ) completely. The values of 
A n and m' n are not used elsewhere in the analysis. 

Our choice of t is made using the "reconstruction" 
scheme [17] [24]: the masses, m n , extracted from the fits to 
the principal correlators, and the Zf extracted from the 
eigenvectors at a single time slice are used to reconstruct 
the correlator matrix using Eq. (16). This reconstructed 



matrix is then compared with the data for t > to, with 
the degree of agreement indicating the acceptability of 
the spectral reconstruction. Adopting too small a value 
of to leads to a poor reconstruction of the data for t > to- 
In general, the reproduction is better as to is increased 
until increased statistical noise prevents further improve- 
ment. The sensitivity of extracted spectral quantities to 
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FIG. 2: The magnitude of matrix elements in a matrix of 
correlation functions, CV,- / ' s/CuCjj , at time-slice 5 is shown 
according to the scale at the lower right. The matrix is for 
the Nucleon H u irrep, with 28 [J = |] operators, 16 [J = |] 
operators and 4 [J = ~] operators. 



the value of to used will be discussed in detail in section 
|VII A[ but in short the energies of low-lying states are 
rather insensitive to to and the reconstruction of the full 
correlator matrix usually is best when to > 7, but not 
too large. 

From the spectral decomposition of the correlator, 
equation [16] it is clear that there should in fact be no 
time dependence in the eigenvectors. Because of states 
at higher energies than can be resolved with dim(C) op- 
erators, there generally is a contribution to energies and 
Z's that decays more rapidly than the lowest mass state 
that contributes to a principal correlator. As for ener- 
gies, we obtain "time-independent" overlap factors, Zf 
from fits of Z"(t), obtained from the eigenvectors, with a 
constant or a constant plus an exponential, in the spirit 
of the perturbative corrections outlined in [14]). 



VI. DETERMINING THE SPIN OF A STATE 

A. Motivation and procedure 

A new method for identifying the spins of excited 
mesonic states was introduced in Refs. [TH] Q~7] . In this 
work, we extend the method to identify spins of excited 
baryonic states. The explanation for the success of the 
new method is that there is an approximate realization 
of rotational invariance at the scale of hadrons in our cor- 
relation functions. There are two reasons for this claim. 
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The first is that there are no dimension-five operators 
made of quark bilinears that respect the symmetries of 
lattice actions based on the Wilson formalism and that do 
not also transform trivially under the continuum group 
of spatial rotations. Thus, rotational symmetry breaking 
terms do not appear until 0(a 2 ) in the action. This argu- 
ment holds even though the action used in this work de- 
scribes an anisotropic lattice. The second reason is that 
our baryon operators are constructed from low momen- 
tum filtered, smeared quark fields, where the smearing is 
designed to filter out fluctuations on small scales. Pos- 
sible divergent mixing with lower dimensional operators 
is then expected to be suppressed. The baryon opera- 
tors then are reasonably smooth on a size scale typical 
of baryons, which is R w 1 fm. With the lattice spacing 
of ~ 0.12 fm, the breaking of rotational symmetry in the 
action or the baryon operators can be small: (a/R) 2 w 
0.015. Of course, these qualitative arguments should be 
backed up by evidence of approximate rotational invari- 
ance in explicit calculation. 



The operators constructed in Section IV using sub- 
duction matrices transform irreducibly under the allowed 
cubic rotations, that is they faithfully respect the sym- 
metries of the lattice. They also carry information about 
the continuum angular momentum, J, from which they 
are subduced. To the extent that approximate rotational 
invariance is realized, we expect that an operator sub- 
duced from spin J to overlap strongly only onto states 
of continuum spin </, and have little overlap with states 
of different continuum spin. In fact this is clearly appar- 
ent even at the level of the correlator matrix as seen in 
Figure [2j Here the correlator matrix for the Nucleon H u 
is observed to be approximately block diagonal when the 
operators are ordered according to the spin from which 
they were subduced. 

To identify the spin of a state, we use the operator 
"overlaps" Zf = (n|O||0) for a given state extracted 
through the variational method presented in the previ- 
ous section. In Figure [3] we show the overlaps for a set 
of low-lying states in the Nucleon H u and Gi u irreps of 
the m,r = 524 MeV 16 3 calculation. The principal corre- 
lators for these states in H u are shown in Figure [T] and 
the corresponding mass spectrum is shown in Figure [9] 
The overlaps for a given state show a clear preference for 
overlap onto only operators of a single spin. While we 
show only a subset of the operators for clarity, the same 
pattern is observed for the full operator set. 

The assignment of spin must hold for states with 
continuum spin J subduced across multiple irreps. In 
the continuum our operators are of definite spin such 
that (0\O J ' M ^J',M') = ZW5 Jt j>5 MtM ; and therefore 
from Eq. [14] the overlap of the subduced operator is 
(a\0 [ $ r \J',M) = Sif^Z^Sj.j,. Only the spin J states 
will contribute, and not any of the other spins present in 
the irrep "A. Inserting a complete set of hadronic states 
between the operators in the correlator and using the 
fact tha t th e subduction coefficients form an orthogonal 

J,M C J,M r r 

= 0A,A'0 r ,r' 



matrix (B7), J2m S A™ S a 



, we thus obtain 



terms in the correlator spectral decomposition propor- 
tional to Z^*Z^ for each A we have subduced into, up 
to discretization uncertainties as described above. Hence, 
for example, a J — | baryon created by a [ J = |] oper- 
ator will have the same Z value in each of the Gi,H, Gi 
irreps. This suggests that we compare the independently 
obtained if-values in each irrep. In Figure|4]we show the 
extracted Z values for negative-parity states suspected of 
being spin | across the H u and Gi u irreps, and of being 
spin | across the Gi u , H u and G2 U irreps. As can be 
seen, there is good agreement of Z values in the different 
irreps that are subduced from spin J, with only small 
deviations from exact equality. 

These results demonstrate that the Z values of care- 
fully constructed subduced operators can be used to iden- 
tify the continuum spin of states extracted in explicit 
computation for the lattices and operators we have used. 

We take the next step and use the identification of the 
components of the spin-J baryon subduced across mul- 
tiple irreps to make a best estimate of the mass of the 
state. The mass values determined from fits to princi- 
pal correlators in each irrep differ slightly due to what 
we assumed to be discretization effects and, in principle, 
avoidable fitting variations (such as the fitting intervals) . 
We follow Ref. [17 and perform a joint fit to the principal 
correlators with the mass being common. This method 
provides a numerical test that the state has been identi- 
fied. We allow a differing second exponential in each prin- 
cipal correlator so that the fit parameters are m n , {m n A } 
and {A^}. These fits are typically very successful with 
correlated \ 2 /Na [ close to 1. suggesting again that the 
descretization effects are small. An example for the case 
of | components identified in Gi u , H u , Gi u is shown in 
Figure [5] When we present our final, spin-assigned spec- 
tra it is the results of such fits that we show. 



B. Additional demonstration of approximate 
rotational invariance 



A further demonstration of approximate rotational in- 
variance is based on the spectrum of energy levels. If 
indeed the mixing between states of different continuum 
spin J is small, then the omission of such coupling should 
not much affect the excited state spectra. That propo- 
sition can be tested by extracting energies using all op- 
erators, and comparing them with the energies obtained 
from only operators subduced from a single J value. If 
approximate rotational invariance were achieved in the 
spectrum, the energies would be nearly the same. As 
an example, we show results for the Nucleon H u irrep, 
in Fig. [6} The left column of the Fig. |6j labelled "all" , 
shows the lowest 12 energy levels obtained from matrices 
of correlation functions using the set of all 48 H u oper- 
ators, and spin identified using the methods previously 
described. The states listed in Figure [3] correspond to a 
few of these 12 levels. The second column shows the low- 
est 6 levels resulting from the variational method when 
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FIG. 3: Histogra ms o f spectral overlaps, Z, are shown for a selection of eight operators (shown at the top using the naming 
convention of Eq. (l3|). The 4 operators labelled with J = | (color red) have subductions only to the H u irrep. The 3 operators 



labelled with J = _ | _ (color green) and the 1 labelled with J = | (color blue) have subductions to both H u and G*2 U irreps. 
Each histogram is labelled by the value of mass m of the state (in units of ma) and has 8 vertical bars showing the relative Z 
values for each of the operators. The data are from the = 524 MeV ensemble and Z's are normalized so that the largest 
value across all states, for a given operator, is equal to 1. The lighter area at the head of each bar represents the one sigma 
statistical uncertainly. Note that for each state only one or two operators have a large relative Z value, and it is the same 
operators appearing across each of the irreps. Note also that nearly the same energies are obtained in H u and G2u irreps for 
state subduced from one J value. 
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FIG. 4: Selected Z values across irreps A„ are shown for states suspected of being J = | (left panel) and | (right panel), based 

on the m w — 524 MeV ensemble. The boxes at the top show the mass for the various states. There are two states of J = | . 

The operators in the left panel, all projected onto J = § , are left to right, Nm ® (\ )m ® ^l}=2 s> ® (I )s ® ^l=2 M' 

® (f + )m ® ^iLi, M > ® (| + )s ® -Dili.M! N ™ ® (|~)m ® -Di=2,S' ® (§~)s ® £>£L 2 ,m- Overlaps of these operators 
after subduction into H u and G2u, agree well for each of the two states shown. The operators in the right panel, all projected 

onto J = |", are left to right, iV M ® (|")m ® D l *l 2M , N M ® (§~)m ® -D£L 2) S> ^ ® (§~)J ®£>ii 2 ,M> ^ ® (|~)s ® ^Lla.M- 
Similarly, the operator overlaps for this state agree well across G\ u , H u , and G2u- 



the operator basis is restricted to only those with contin- 
uum spin J = |. Similarly, the third shows the lowest 4 
levels when the operator basis is restricted to only those 
with continuum spin J = |, and the last column are 
the lowest two levels when the basis is restricted to only 
,/=| operators. 



The results are striking. We see that the masses of the 
levels in each of the restricted bases agree quite well with 
the results found in the full basis. The agreement is quite 
remarkable because one expects operators in the H u irrep 
that can couple with the ground state (the lowest J = I 
state), will allow for the rapid decay of correlators down 
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FIG. 5: Fit to the three subduced principal correlators of lowest lying | Nucleon using a common mass. Results are from 
the m,r = 524 MeV ensemble. These levels correspond to the sixth excited Gi u , the seventh H u and the fourth Gi u . Plotted 
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FIG. 6: Extracted Nucleon H u mass spectrum for various 
operator bases, (a) Full basis (dim = 48), (b) Only J = | 
operators (dim = 28), (c) Only J = ~ operators (dim = 16), 
(d) Only J = ~ operators (dim = 4). Results are from the 
ra,r = 524 MeV ensemble. The dimensionality of the operator 
basis in each J is shown in Table [Alj] 
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FIG. 7: Extracted Nucleon H u mass spectrum as a function 
of to- Horizontal bands to guide the eye. For clarity of pre- 
sentation, only one of the highly excited J = | and J = | 
levels are shown. 



analysis, and the number of distillation vectors. We will 
use the Nucleon H u and the Delta H g irreps in the m, = 
524 MeV dataset to demonstrate our findings. 



to the ground state as a function of time, t. However, the 
higher-energy spin | and \ states do not show such a de- 
cay; we obtain good plateaus in plots like those in Fig. [T] 
These results provide rather striking demonstration for 
the lack of significant rotational symmetry breaking in 
the spectrum. 



VII. STABILITY OF SPECTRUM 
EXTRACTION 

In this section we consider to what extent the extracted 
spectrum changes as we vary details of the calculation, 
such as the metric timeslice, io; used in the variational 



A. Variational analysis and to 

Our fitting methodology was described in Section |V| 
where reconstruction of the correlator was used to guide 
us to an appropriate value of to. As seen in Figure [7J for 

^ 7, the low-lying mass spectrum is quite stable with 
respect to variations of to . This appears to be mostly due 
to the inclusion of a second exponential term in Eq. ( 17 1, 



which is able to absorb much of the effect of states out- 
side the diagonalization space. The contribution of this 
second exponential typically falls rapidly with increasing 
to both by having a smaller A and a larger m' . 

Overlaps, Z l n = (n|O]|0>, can show more of a sensitiv- 
ity to to values being too low, as was found also in the 
analysis of mesons, Ref. [I?] . 
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FIG. 8: Extracted the Delta H g mass spectrum as a func- 
tion of number of distillation vectors in the m w = 524 MeV 
dataset. For N < 32, spin identification is lost for the highest 
J = | + and J = | + levels. 



angular dependence of the wavefunction over the typical 
size of a hadron. The results for baryons presented here 
are consistent with these observations. 

In summary one is limited as to how few distillation 
vectors can be used if one requires reliable extraction 
of high-spin states. The results shown here suggest 32 
distillation vectors on a 16 3 lattice is the minimum, so at 
least 64 distillation vectors on a 20 3 lattice are likely to 
be required. 



VIII. RESULTS 

Our results are obtained on 16 3 x 128 lattices with 
pion masses between 396 and 524 MeV. More complete 
details of the number of configurations, time sources and 
distillation vectors used are given in Table |l} T he full 
basis of operators in each irrep listed in Table [Hi] is used 
for the variational method construction. 



In summary it appears that variational fitting is 
reliable provided to is "large enough". Using two- 
exponential fits in principal correlators we observe rel- 
atively small to dependence of masses, but more signif- 
icant dependence for the Z values which we require for 
spin-identification. 



B. Number of distillation vectors 

The results presented so far are based on the analysis 
of correlators computed on 16 3 lattices using 56 distilla- 
tion vectors. We consider how the determination of the 
spectrum varies if one reduces the number of distillation 
vectors and thus reduces the computational cost of the 
calculation. This is particularly important given that, 
as shown in |19j . to get the same smearing operator on 
larger volumes one must scale up the number of distil- 
lation vectors by a factor equal to the ratio of spatial 
volumes. To scale up to a 32 3 lattice this would require 
56 x (f|) 3 = 448 vectors which is not currently a realiz- 
able number without using stochastic estimation [25j . 

In figure [8] we show the low- lying part of the extracted 
Delta H g spectrum on the m w = 524 MeV lattice as 
a function of the number of distillation vectors used in 
the correlator construction. It is clear that the spectrum 
is reasonably stable for TV > 32 but that the spectrum 
quality degrades rapidly for fewer vectors. In particular, 
the methods for spin identification fail for the highest 
J p = | + state as well as the first J p = ^ + state. 

The need for a large numbers of distillation vectors 
has been discussed in Ref. [17] for the case of isovector 
mesons. The conclusion drawn is that to describe high 
spin hadrons having large orbital angular momentum, 
one needs to include sufficient vectors to sample the rapid 



A. m v = 524 MeV results 

Using the variational method outlined in Section [VJ 
the spectrum of energies in each lattice irrep for the Nu- 
cleon and for the Delta are shown in Figures [9] and [TUJ re- 
spectively. By applying the spin-identification procedure 
described above, we can subsequently associate contin- 
uum J p labels to each of these states, as labelled on the 
figures. For the remainder of this paper, we will there- 
fore label the energies by their assigned continuum spins; 
these are shown for the lattices at m, = 524 MeV in 
Figure [IT] with the energies obtained from a joint fit to 
the principal correlators as described in Section |Vj 
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FIG. 9: Extracted Nucleon spectra by irrep for m w = 
524 MeV. Colors are black (J = |), red (§), green (§), 
blue (|). Masses are shown in ratios of the Q baryon mass. 

There are several notable features in these spectra. 
As we will discuss, the patterns of states have a good 
correspondence with single-hadron states as classified by 
SU(6)<E>0(3) symmetry. The numbers of low-lying states 
in each J p are similar to the numbers obtained in the 
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FIG. 10: Extracted Delta spectra by irrep for = 524 MeV. 
Colors are black ( J = |), red (|), green (|), blue (|). Masses 
are shown in ratios of the Q, baryon mass. 



non-relativistic quark model which is a particular real- 
isation of the symmetry above (e.g., &7\). For the 
purposes of these comparisons, it is helpful to introduce 
a spectroscopic notation: X 2S+1 L 7r J p , where X is the 
Nucleon N or the Delta A, S is the Dirac spin, L = S, 
P, D,. . . denotes the combined angular momentum of the 
derivatives, 7r = S, M, or A is the permutational symme- 
try of the derivative, and J p is the total angular momen- 
tum and parity. This notation also is used in Table |TV} 
which we discuss now. 

In the negative parity N* spectrum, there is a pat- 
tern of five low- lying levels, consisting of two lev- 
els, two iV| levels, and one V| level. The triplet 
of higher levels in this group of five is nearly degen- 
erate with a pair of A A and A| levels. This pat- 
tern of Nucleon and Delta levels is consistent with an 
L = l~ P-wave spatial structure with mixed symme- 
try, Pm . As shown in Table |IV[ the same numbers of 
states are obtained in the SU(6) ® 0(3) classification for 
the negative-parity Nucleon and Delta states constructed 
from the "non-relativistic" Pauli spinors as we find in the 
lattice spectra. The lowest two N*~ states are dominated 
by operators constructed in the notation of Eq. [13] as 
N M (g)(S = ^ + )m<E)(L = 1 _ ) M ->• J p = |~ and §~, while 
the three higher N*~ levels are dominated by operators 
constructed according to A^m ® (S = | + )s ® (L = 1 _ )m 

Similarly, the low-lying Delta 



with J p = A 



I andf 



levels are consistent with a A 1 , and A| assignment. 
There are no low-lying negative-parity S = § Delta states 
since a totally symmetric state (up to antisymmetry in 
color) cannot be formed. Consequently, there is no low- 
lying A | , which agrees with the lattice spectrum. In 
the non-relativistic quark model [26] . a hyperfine contact 
term is introduced to split the doublet and quartet states 
up and down, respectively, compared to unperturbed lev- 
els and the tensor part of the interaction provides some 
additional splitting. The result is that the doublet Delta 
states are nearly degenerate with the quartet Nucleon 
states as is observed in the lattice spectra, Fig. 11 In 
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TABLE IV: Subset of the operator basis classified by SU(6)® 
0(3) multiplets and total spin and parity J p . The entries 
are the orbital angular momentum structures outlined in Sec- 



tion 



IV 



and Appendix [X] that contribute within each J p . The 
operators listed here are all from the first embedding in Dirac 
spin in Table |VI| and correspond to only upper components 
(referred to as "non-relativistic"). A spectroscopic notation 
of X 2S+1 L 7T J P is used, where X = N or A, S is the Dirac 
spin, L = S, P, D,... is the combined angular momentum of 
the derivatives, n = S, M, or A is the permutational symme- 
try of the derivatives, and J p is the total angular momen- 
tum and parity. In Sections |IVj|VI| an operator notation like 



(Am ® ( 



3+\ 1 
2 JS 



>D ll] 



was used which in this spectro- 



scopic notation would be N Pm§ • Dimensions and parities 
of SU(6) <g> 0(3) representations are listed in column 4 for 
the doublet spin states of column 2, and in column 5 for the 
quartet spin states of column 3. The number, n, of opera- 
tors for each J p is listed in the rightmost column. This same 
spectroscopic notation and classification of spatial structure 
is also used for comparisons with models where L represents 
the orbital angular momentum. 



the language of SU(6) ® 0(3), these low-lying N and A 
states constitute the strangeness zero part of a [70, 1 _ ] 



multiplet, as indicated in Table IV 

In the positive parity sector, there also are inter- 
pretable patterns of lattice states in the range m/rriQ ~ 

1 + ^ H- 

1.3 — 1.5. There are four N ^ levels, five A~| levels, 
three 7V| levels and one N-^ level, which are the same 
numbers of levels for each J p as in Table llVI The lattice 
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FIG. 11: Spin-identified spectrum of Nucleons and Deltas from the lattices at = 524 MeV, in units of the calculated f2 
mass. 
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FIG. 12: Spin-identified spectrum of Nucleons and Deltas from the lattices at m n = 396 MeV, in units of the calculated Q 
mass. 



spectra also have two A^ levels, three A| levels, two 
A| + levels, and one A| + level, which are the same num- 
bers of levels for each J p as in Tabic llVl In this case we 
are considering the multiplets [70,0+], [56,2+], [70,2+], 
[20, 1 + ] and a radially excited [56, + ] and within non- 
relativistic qqq constituent quark models |27] , the mass 
eigenstates are admixtures of these basis states. 

In general for both the Nucleon and Delta spectrum, 
there are reasonably well-separated bands of levels across 
the range of J values, alternating in parity, with each 
band higher in energy than the previous one. We remark 
that there are no obvious patterns of degenerate levels 
with opposite parities for the same total spin, J as in 
Ref. M- 



The discussion up to now has focused on observables; 
namely, the level energies. More information about the 
internal structure of the lattice states can be obtained by 
analyzing the spectral overlaps. We remind the reader 
that the full basis of operators listed in Table |III| is 
used within the variational method. We find that only a 
few operators have large overlaps in the lowest negative- 
parity Nucleon levels of the G lu , H u and G^u irreps, 
and these are the subduced versions of "non-relativistic" 
operators. From the construction of operators in Ap- 
pendix |Aj those featuring only upper components in spin 
are in the first embedding of Dirac spin in Table |VI| 
There is a one-to-one correspondence between these op- 



erators, and those listed in Table IV and we will adopt 
the spectroscopic names as a shorthand. This spectro- 
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FIG. 13: Spectral overlaps for the negative parity Nucleon 
operators in Table |IV| within the lowest lying states for the 
m,r = 524 MeV ensemble. The operators are subduced into 
the Gi u , H u and Gi u irreps, and their spectral overlaps are 
shown for the lowest lying states with the masses listed in 
units of the fl baryon mass. In the spectroscopic notation, 
these are the N 2 Pm s 



and 7V 4 Pm| subduced into Gi u , the 



JV 2 P M § andAr 4 P M -i 
duced into both H u and G2u- 
one operator is dominant. 



g sub- 
As can be seen, for each level 



scopic notation is used in Table |TVj to identify operators, 
but is also applicable for identification of states. (De- 
tails of the derivative construction for the operators can 
be found in Appendix |A 3| ) The Nucleon (and Delta) 
J = \ i f and | operators all feature one derivative 
in P M coupled to either a spin S = | (doublet) or S = | 
(quartet). Spectral overlaps (Z) for these operators can 
be directly compared since they all have a consistent nor- 
malization. An S = I spin coupled to one derivative can 



be projected to either J 



J 



or J 



The first 



and second constructions are subduced into the G\ u and 
H u irreps, while the J = § construction is subduced into 
the H u and G^u irreps. 

In Figure [13} we compare these non-relativistic opera- 
tor overlaps for the lowest lying states and across irreps. 
We see that one operator is dominant for each state, with 
magnitudes that are roughly consistent across all irreps. 
In particular, the constructions coupled to different J 
have similar magnitudes. These results suggest that these 



low-lying levels form a multiplet with little mixing among 
the states. In the language of SU(6) (g> 0(3) (spin-flavor 
and space), these states, and their negative parity Delta 
partners, are part of a [70, 1 _ ] multiplet. 

Similarly, we can examine the first group of excited 
positive-parity levels in the Nucleon channel that cluster 
around m/mn ~ 1.4. Again, we find that the positive 
parity Nucleon operators from Table |IV| are dominant 
in the spectral overlaps. One of them has a quasi-local 
structure (indicated as Ss for no derivatives) with J = h, 
while the others involve two derivatives coupled either 
to L = 0, 1, or 2, and are labeled as 5m, Pa or D$ 
and Dm- We find that there is not a unique mapping 
of each state to one particular operator, rather each op- 
erator contributes in varying magnitude to each state, 
indicating significant mixing in this basis. 

These results, along with the observation that the 
numbers of states are consistent with the numbers of 
non-relativistic operators in each J + , suggest that this 
band of positive-parity states belongs to more than one 
SU(6)(&0(3) multiplet, with now mixing among the mul- 
tiplets mentioned before; namely, the [70, + ], [56, 2 + ], 
[70, 2+] [20, 1+] and a radially excited [56, 0+]. It is no- 
table that there is overlap with all the allowed L^J P 
multiplets with L < 2, and in particular, there is mixing 
with the [20, 1 + ] multiplet. There does not appear to be 
any "freezing" of degrees of freedom as suggested in some 
diquark models (for some reviews see Refs. [2S1I30])- We 
will return to this point in the summary. 

As we move up to the second excited negative-parity 
band in the Nucleon and Delta channels, we find the 
"non-relativistic" Pm operators discussed previously do 
not feature prominently in the spectral overlaps in these 
higher lying excited states. Instead, operators using two 
derivatives together with one quark having a lower com- 
ponent (ft = — ) Dirac structure appear, some of which 
appear in Figures [3] and [4] The lower component spinor 
contributes a factor — 1 to the parity. Similarly, those op- 
erators that featured prominently in the first excited pos- 
itive parity band do not appear significantly in the second 
excited positive band in the Nucleon and Delta chan- 
nels. Instead, operators involving two derivatives and 
two lower-component Dirac spinors appear. The lower 
components effectively bring in angular momentum, but 
are not equivalent in spatial structure. To adequately 
resolve the internal structure of these higher lying states 
will require the introduction of operators featuring three 
and four derivatives. 



B. Quark mass dependence 

At the two lighter pion masses, = 444 and 396 
MeV, we find the spectra for Nucleons and Deltas, clas- 
sified by irreps, to be qualitatively similar to the heaviest 
pion mass case. The spin identification techniques de- 
scribed above are successful, and the resulting spectrum, 



identified by spin is shown in Fig. 12 for the ensemble at 
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FIG. 14: Lightest few Nucleon J 
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states. Also shown in 7V| + is the threshold for Ntt and Ntttt. 
The influence in the spectrum from these thresholds is com- 
plicated by the use of a finite spatial cubic box. Further dis- 
cussion is in Section IPXl 
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FIG. 15: Lightest few Nucleon J = 
states. Also shown in 7V| is the threshold for Ntt and Ntttt. 
The influence in the spectrum from these thresholds is com- 
plicated by the use of a finite spatial cubic box. Further dis- 
cussion is m Section ILXl 



TOtt = 396 MeV, normalized in units of the measured f2 
baryon mass determined at this quark mass. 

The Nucleon and Delta spectrum is qualitatively sim- 
ilar to that determined at the heaviest pion mass - the 
TOtt = 524 MeV lattices - albeit typically at smaller mass 
in units of the Q baryon mass. There are the same num- 
ber and pattern of low lying negative parity Nucleon and 
Delta states, as well as in the first excited band of positive 
parity states. Again, the low lying negative parity Delta 
states are slightly higher than the corresponding Nucleon 
states. There is a slightly enhanced splitting for the first 
excited band of positive parity Delta states across J + , as 
well as more mixing among the positive parity Nucleon 
states around m/mn ~ 1.3. 

In Figures [14] - [TTJ we show extracted state masses as 
a function of la = § fe^4i which we use as a proxy for 
the quark mass [11] for the three quark masses listed in 
Table |l] The state masses are presented via ^^ m Q hyS ' ■ 
The ratio of the state mass (ma) to the fi-baryon mass 
computed on the same lattice removes the explicit scale 
dependence and multiplying by the physical Sl-baryon 
mass conveniently expresses the result in MeV units. 
This is clearly not a unique scale-setting prescription, but 



it serves to display the data in a relatively straightfor- 
ward way. We remind the reader that the data between 
different quark masses are uncorrelated since they follow 
from computations on independently generated dynami- 
cal gauge- fields. 

In Figure [14] we show the mass of some of the lowest 
identified positive parity levels. In Figure [15] we show 
the lowest few negative parity Nucleon levels. To help 
resolve near degeneracies, we slightly shift symbols at 
the same pion mass horizontally in Iq. In some cases, 
for comparison, we plot the mass of the lightest A?V and 
Ntttt thresholds - the mass follows from the simple sum 
of the extracted masses on these lattices. In general, the 
lattice levels decrease with the quark mass. There is no 
observed dramatic behavior, for example, from crossing 
of thresholds. 

Notable in these plots is the clustering of bands of lev- 



els as seen in Figures 11 and 12 and described above. 
In the first excited positive parity Nucleon band, there 
is a tendency for the levels across J + to cluster, and is 
most readily apparent in the J = ^ + to cluster around 
m/rriQ ~ 1.3 as the pion mass decreases. Another no- 
table feature in these figures is the absence of an excited 



Nucleon J 



i + 



state comparable, or slightly below, the 
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FIG. 16: Lightest few Delta J 



f + , § + , | + , and | + states. 
Also shown in A| + is the threshold for Ntt and Nitty. The 
influence in the spectrum from these thresholds is complicated 
by the use of a finite spatial cubic box. Further discussion is 
in Section ITXl 



FIG. 17: Lightest few Delta J = 



| , | , | , and ~ states. 
Also shown in A | is the threshold for Ntt and Ntttt. The 
influence in the spectrum from these thresholds is complicated 
by the use of a finite spatial cubic box. Further discussion is 
in Section llXl 



lowest lying ~ level as reported in the PDG [31] • Simi- 
larly, there is no low lying Delta J = | state comparable 

in mass to the | . The ordering of the low-lying states 
in the Nucleon spectrum, and in particular the spectrum 
of the channel, has been the subject of much effort 

in lattice QCD. 



C. Comparisons 



In Figure 18 we show a comparison of our result for 
the Nucleon J = \ + spectrum with other calculations 
in full QCD. Those shown in grey are from Ref. [32] 
using 2 © 1 dynamical quark flavors, while those in or- 
ange are from Ref. jl] using two dynamical light-quark 
flavors. Each of these calculations employs a different 
means to set the lattice scale, and we have made no at- 
tempt to resolve these calculations to a common scale- 
setting scheme. Rather our aim in this discussion is to 
compare the pattern of states that has emerged in each 
calculation, and to provide a possible explanation for the 
differences. 

References pS] and |3] find only two excited levels be- 



low 2.8 GeV, notably isolated from one another in the 
case of Ref. [32], at pion masses comparable to those 
in this study. At the lightest quark mass reported in 
this work, there are jour nearly degenerate excited states 
found at approximately 2.2 GeV, and three nearly de- 
generate states near 2.8 GeV. A possible explanation for 
the discrepancy in the number of levels is the operator 
constructions used. Refs. [H |32] use a basis of local or 
quasi-local operators, without, for example, derivatives, 
but with multiple smearing radii. These operators, which 
have an S'-wave spatial structure concentrated at the ori- 
gin, can have overlap with radial excitations of S'-waves, 
but will have limited overlap with higher orbital waves. 
The results presented here suggest the observed excited 
J ~ \ + states are admixtures of radial excitations as well 
as Z?-wave and anti-symmetric P-waves structures, and 
the inclusion of operators featuring such structures is es- 
sential to resolve the degeneracy of states. The impact of 
an incomplete basis of operators will be addressed more 
in Section HXl 

The authors of Ref. [32] note the large drop in the ener- 
gies of the first and second excited states at their lightest 
pion mass, and ascribe this to the emergence of a light 
"Roper" in their calculation. While we do not reach cor- 
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FIG. 18: Comparison of results for the Nucleon J = | + channel. The results shown in grey are from Ref. [32], while those in 
orange are from Ref. [4[. Note that data are plotted using the scale-setting scheme in the respective papers. Results from this 
paper are shown in red (the ground state), green and blue. At the lightest pion mass, there is a clustering of four states as 
indicated near 2 GeV, while there are three nearly degenerate states 2.7 GeV. Operators featuring the derivative constructions 
discussed in this paper feature prominently in these excited states. 



respondingly light pion masses and large volumes in this 
study, the work presented here clearly shows the need for 
a sufficiently complete basis of operators before a faithful 
description of the spectrum can emerge, and the identi- 
fication of the Roper resonance be warranted. Indeed, 
even our work remains incomplete since we have entered 
a regime of open decay channels, discussed next, requir- 
ing the addition of multi-particle operators to the basis. 



IX. MULTI-PARTICLE STATES 

In the previous section we presented the extracted 
spectra from calculations with three different light quark 
masses. In each case we were able, using the operator 
overlaps, to match states across irreps that we believe 
are subduced from the same continuum spin state. This 
suggests an interpretation of the spectrum in terms of 
single-hadron states, while in principle our correlators 
should receive contributions from all eigenstates of finite- 
volume QCD having the appropriate quantum numbers. 
This includes baryon-meson states which in finite volume 
have a discrete spectrum. Where are these multi-particle 
states? 

This issue was investigated in the meson sector [T7] 
where the spectrum was compared between multiple vol- 
umes and on multiple mass data sets. In particular, un- 
der change of volume, the extracted spectrum did not 
resemble the changing pattern of levels one would expect 
from two-meson states, but rather was largely volume- 



independent. As such, the interpretation of the observed 
levels was that of a single particle spectrum. In this work, 
some initial investigations were made with a 20 3 x 128 lat- 
tice at the lightest pion mass available, and similar obser- 
vations are made; namely, the spectrum between different 
volumes does not change substantially. We are thus led 
to interpret the spectrum is terms of single hadron states. 
The subsequent observations we make are quite similar 
to those made for mesons [17] . 

The overlap of a localized three-quark operator onto a 
baryon-meson state will be suppressed by l/\/V, where 
V is the lattice volume, if the operator creates a reso- 
nance with a finite width in the infinite volume limit. 
This fall-off is matched by a growth in the density of 
states with the volume and the resonant state thus main- 
tains a finite width as the mixing with each discrete state 
falls. The simulations in this study are carried out in cu- 
bic volumes with side-lengths ~ 2 fm, which might be 
sufficiently large that the mixing between one of the low- 
lying two-particle states and a resonance is suppressed 
sufficiently for it to be undetectable with the three-quark 
operator basis. 

Even if the mixing between localized single-hadron 
states and baryon-meson states to form resonance-like fi- 
nite volume eigenstates is not small, there still remains a 
practical difficulty associated with using only three-quark 
operators. In this case the state can be produced at the 
source time-slice through its localized single-hadron com- 
ponent, while the correlator time dependence obtained 
from e~ Ht will indicate the mass of the resonant eigen- 
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state. Consider a hypothetical situation in which a sin- 
gle baryon-meson state, denoted by |2), mixes arbitrarily 
strongly with a single localized single-hadron state, |1), 
with all other states being sufficiently distant in energy 
as to be negligible. There will be two eigenstates 

| a) = cos0|l) +sin6»|2) 
\b) = -sin0|l) +cos(9|2), 

with masses m ai m\,. At the source (and sink) only the 
localized single-hadron component of each state overlaps 
with the operators in our basis and hence the overlaps, 
Zp b = (ft, b|O]|0), will differ only by an overall multi- 
plicative constant, Zf = cos 6»zf \ Z\ = -sinflzj 1 *. As 
such the eigenvectors v a ,v b point in the same direction 
and cannot be made orthogonal. Thus the time depen- 
dence of both states will appear in the same principal 
correlator as 

A(t) ~ A a e- m « {t - to) + A b e- mb(t - to) + ... 

Since m a and m b most likely do not differ significantly 
(on the scale of a^ 1 ) it will prove very difficult to ex- 
tract a clear signal of two-exponential behavior from the 
principal correlator. This is precisely why the varia- 
tional method's orthogonality condition on near degen- 
erate states is so useful, but we see that it cannot work 
here and we are left trying to extract two nearby states 
from a \ 2 St to time-dependence. Typically this is 
not possible and reasonable looking fits to data are ob- 
tained with just one low-mass exponential. This is anal- 
ogous to the interpretation made in the comparison be- 
tween our computed spectrum, and those computed 
using only rotationally-symmetric smeared (quasi-)local 
sources [H |32], shown in Figure 18 our results showed 
a cluster of four near-degenerate states, while the other 
analyses showed one or perhaps two since all four states 
would only couple through their S-wave components. 

Thus, in some portions of the extracted spectrum, we 
might be observing admixtures of "single-particle" and 
baryon-meson states. A conservative interpretation then 
of our spectrum is that the mass values are only accurate 
up to the hadronic width of the state extracted, since 
this width is correlated with mixing with baryon-meson 
states via a scattering phase-shift. 

In order to truly compare with the experimental situa- 
tion, we would like to explicitly observe resonant behav- 
ior, thus to obtain a significant overlap with multi-hadron 
states, we should include operators with a larger number 
of fermion fields into our basis, and in particular, multi- 
particle operators. The construction of single-meson and 
single baryon operators of definite continuum helicity and 
subduced into the 'in flight' little-group irreps can be 
done using the tables in [33J GH] • Spin identification is 
possible and will be reported in future work. These in- 
flight operators can be used in two-particle constructions, 
and along with single particle operators, provide a far 
more complete determination of the excited levels in an 
irrep. 



As shown by Luscher[35], one can map these discrete 
energy levels onto the the continuum energy dependent 
phase shift within a partial wave expansion, including the 
phase shift for higher partial waves. Such a technique 
was recently used to determine the L — and L = 2 
phase shift for non- resonant 1 = 2 tttt scattering [36j . 
The mapping of the phase shift is both volume and irrep 
dependent. This is the origin of the cautionary remarks 
in the captions of Figures 14 - fTTJ Namely, the location 



of the threshold energies are in fact irrep dependent and 
not solely determined by the energy of the continuum 
states. 

With suitable understanding of the discrete energy 
spectrum of the system, the Luscher formalism can be 
used to extract the energy dependent phase shift for a res- 
onant system, such as has been performed for the 1=1 
p system [37] . The energy of the resonant state is deter- 
mined from the energy dependence of the phase shift. It 
is this resonant energy that is suitable for chiral extrap- 
olations. 

Annihilation dynamics feature prominently in resonant 
systems, and these dynamics arise from quark discon- 
nected diagrams in multi-particle constructions. Utiliz- 
ing the techniques developed recently for the study of 
isoscalar systems [35], distillation, and stochastic vari- 
ants [25], can be used for the efficient numerical eval- 
uation of multi-particle systems. 



X. SUMMARY 

We have described in detail our method for extracting 
a large number of Nucleon and Delta excited states using 
the variational method on dynamical anisotropic lattices. 
Key to the success of the method has been the use of 
a large basis of carefully constructed operators, namely 
all three quark baryon operators consistent with classical 
continuum symmetries, and with up to two derivatives, 
that are subsequently projected onto the irreducible rep- 
resentations of the cubic group. We have exploited the 
observed approximate realization of rotational symme- 
try to devise a method of spin identification based on 
operator overlaps, enabling us to confidently assign con- 
tinuum spin quantum numbers to many states. We have 
demonstrated the importance of having a suitable oper- 
ator basis with overlap to all the continuum spin states 
that contribute to the spectrum. We have then demon- 
strated the stability of the spectra with respect to chang- 
ing the number of distillation vectors and the details of 
the variational analysis. We have successfully applied 
these techniques at one lattice volume with three differ- 
ent light quark masses. We are able to reliably extract 
a large number of excited states with J p ranging from 
J = \ up through J = | in both positive and negative 
parity. These are the first lattice calculations to achieve 
such a resolution of states in the baryon sector with spin 
assignments and J > § . 

We find a high multiplicity of levels spanning across J p 
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which is consistent with SU(6) <8 0(3) multiplet count- 
ing, and hence with that of the non-relativistic qqq con- 
stituent quark model. In particular, the counting of lev- 
els in the low lying negative parity sectors are consistent 
with the non-relativistic quark model and with the ob- 
served experimental states [31] . The spectrum observed 
in the first excited positive parity sector is also consis- 
tent in counting with the quark model, but the compar- 
ison with experiment is less clear with the quark model 
predicting more states than are observed experimentally, 
spurring phenomenological investigations to explain the 
discrepancies (e.g., see Refs. [271l29TJ3"Ill3l)H4"I]). 

We find that each of the operators in our basis fea- 
tures prominently in some energy level, and there is sig- 
nificant mixing among each of the allowed multiplets, in- 
cluding the 20-plet that is present in the non-relativistic 
qqq quark model, but does not appear in quark-diquark 
models 29 , and in particular Ref. [42]. This adds fur- 
ther credence to the assertion that there is no "freezing" 
of degrees of freedom with respect to those of the non- 
relativistic quark model. These qualitative features of 
the calculated spectrum extend across all three of our 
quark-mass ensembles. Furthermore, we see no evidence 
for the emergence of parity-doubling in the spectrum [28]. 

We have argued that the extracted spectrum can be 
interpreted in terms of single-hadron states, and based 
on investigations in the meson sector [T7] and initial in- 
vestigations of the baryon sector at a larger volume, we 
find little evidence for multi-hadron states. To study 
multi-particle states, and hence the resonant nature of 
excited states, we need to construct operators with a 
larger number of fermion fields. Such constructions are 
in progress, and we believe that the addition of these op- 
erators will lead to a denser spectrum of states which can 



be interpreted in terms of resonances via techniques like 
Liischer's and its inelastic extensions [43 . 

The extraction and identification of a highly excited, 
spin identified single-hadron spectrum, represents an im- 
portant step towards a determination of the excited 
baryon spectrum. The calculation of the single-baryon 
spectrum including strange quarks is ongoing. Combin- 
ing the methods developed in this paper with finite vol- 
ume techniques for the extraction of phase shifts, future 
work will focus on the determination of hadronic reso- 
nances within QCD. 
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Appendix A: Construction of flavor /spin/space 
symmetric states 

Consider the construction of sets of definite symme- 
try for three objects that can be labeled by x, y, and 
z, where the first object is labeled by x, the second is 
labeled by y and the third is labeled by z. There are in 
general, four definite symmetry combinations: symmet- 
ric, mixed-symmetric, mixed anti-symmetric, and totally 
antisymmetric, denoted by S, MS, MA, A. 



Let symmetry projection operator S-z 



be 



x y z 

K %' y' * J 

defined so that its action on a generic object with la- 
bels x' , y' and z' is to create a superposition of objects, 
denoted by {xyz}-£, with symmetry Z of their labels as 
follows, 



{xyz} T = ^2 S T 

x'y'z' 




{x'y'z'}. 



(Al) 



Permutation operators St can be inferred from their 
action on an object with three labels to produce the four 
allowed symmetry combinations as follows, 

{xyz} s = N s [{xyz} + {yxz} + {zyx} + {yzx} 
+ {xzy} + {zxy}], 
{xyz} M s = N MS [{xyz} + {yxz} + {zyx} + {yzx} 

- 2{xzy} - 2{zxy}], 
{xyz} M A = N M a [{xyz} - {yxz} + {zyx} - {yzx}] , 

{xyz} A = N A [ - {xyz} + {xzy} - {yzx} + {yxz} 

- {zxy} + {zyx}] . (A2) 

If two or more labels are the same, then equivalent terms 
must be combined and the normalization constants ad- 
justed to give an appropriate normalization, e.g., for or- 
thonormal quantum states N$ = when x, y and z 
all are different. By convention, mixed symmetries are 
MS or MA according to whether the first two labels are 
symmetric or antisymmetric. 

The use of projection operators allows the same con- 
structions to be applied to the labels of operators as are 
applied to the labels of the states created by the oper- 
ators. They also provide the basis for a straightforward 
computational algorithm that yields the desired super- 
positions. 

Baryon operators have sets of labels for flavor, spin, 
and spatial arguments that transform independently, 
therefore as direct products. Each of these sets of la- 
bels can be arranged according to Eq. (A2). Then the 



symmetries of the sets of flavor labels must be combined 
with the symmetries of the sets of spin and spatial labels 
in order to make an overall symmetric object as discussed 
in the text. The general rules for combining direct prod- 
ucts of objects with independent labels sets 1 and 2 to 
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make overall symmetries of the combined sets of all la- 
bels, denoted by 1,2, are as follows, 

{l}s{2} s = {1, 2} s , {1}s{2}ms = {1, 2}ms, 
{1}s{2}ma = {1, 2}ma, {1}s{2}a = {1, 2} A , 

{1} A {2} S = {1, 2} A , {1} A {2} MS = {1, 2}ma, 
-{1}a{2} ma = {1, 2} MS , {1}a{2}a - {1, 2} s , 



71' 



{1}ms{2}ms + {1}ma{2}ma) 
{1}ms{2}ms + {1}ma{2}ma) 
{1}ms{2}ma + {1}ma{2}ms) 
{1}ms{2}ma + {1}ma{2}ms) 

1. Dirac spin 



{1,2} S , 
{!■ 2}ms, 
{1, 2}ma, 
{1,2}a 



(A3) 



The construction of Dirac spin representations follows 
Table 14 in Appendix B of Ref. [7] . A Dirac spinor with 
an index that takes four values is formed from direct 
products of two spinors that have indices that take two 
values: one for ordinary spin (s-spin) and the other for 
intrinsic parity (p-spin) . The s-spin indices are s = + for 
spin up and s = — for spin down while the p-spin indices 
arc p = + for positive intrinsic parity and p = — for neg- 
ative intrinsic parity. These p- and s-spin indices deter- 
mine the Dirac spin index (in the rest frame) as shown 
in Table |Vj based on the Dirac-Pauli representation of 
Dirac matices. 

There are eight s-spin states for a three-quark baryon 
and they are obtained by projecting an arbitrary state to 
combinations with good symmetry. The resulting states 
follow from Eq. (A2| with x, y and z taking two values, 
+ and — , which produces four symmetric (S) states with 
total spin |, 



I 2 ' 

I 2 ' +2~)s 

I 2 ' 

I 3- 
I 2 ' 



|3 ,+|> s = |+ + + ) s , 

! ; -i++-> s . 

2~)s = H )s' 

■§>s = !--->s> ( A4 ) 
two mixed-symmetric (MS) states with total spin |, 

I 2> +2~)mS = ~H ~i ^ )mS' 

If' - 2~)ms = — I ' )ms' 
and two mixed-antisymmetric (MA) states with total spin 



1 + 1 ) 

2> ' 2/MA 

I _I\ 
2' 2/MA 



+ - + ) 



Dirac 


P s 


1 


+ + 


2 


+ - 


3 


- + 


4 





TABLE V: Mapping of Dirac spin indices to p and s labels. 



The states on the right side are normalized states of def- 
inite symmetry. The eight p-spin states take exactly the 
same form as the s-spin states. 

The 64 Dirac spin labels of three quarks are obtained 
from direct products of p-spin and s-spin states of the 
quarks. The possible symmetries of the Dirac spinors 
arc obtained from the multiplication rules in Eq. |A3| to- 
gether with Table [V] which shows how each quark's Dirac 
index is determined by its s-spin and p-spin indices. Ex- 
amples of the construction are given in Ref. [7]. Equal 
numbers of positive- and negative-parity states are al- 
ways produced. The octahedral irrep of the product is 
G\ for s-spin \ and H for s-spin | representations of 
SU(2). Table [Vl] shows the irreps of St/ (2) that are pro- 



duced. The number of embeddings of each irrep also is 
listed. Parities of the states are determined by products 
of the p-spins of the three quarks. 



2. Flavor 

The flavor states with wel l-defined symmetries also are 
constructed using Eq. (A2). For the purposes of dis- 



MA' 
MA' 



(A6) 



cussion, we will consider the flavor symmetric limit, i.e., 
degenerate quark masses. There are three separate rep- 
resentations - the symmetric (S) states are the SU(3) 
decuplet, 1 10) s , the mixed-symmetric (MS, MA) states 
are the £77(3) octets, |8) MSMA , an d the antisymmetric 
(A) state is the 5/7(3) flavor-singlet state, |l) A - The 
A + is \udu)s, while the proton is \ucIu}ma,ms- The E° 
is |mc?s)ma,ms- The octet, singlet and decuplet construc- 
tions are shown in Table IVIII 

Within the SU (3) flavor representations, we also have 
SU(2) isospin states. In the construction that follows, 
it is straightforward to generalize to the case of broken 
SU(3). The combinations of symmetry states remain 
valid, however, there are some new states, such as E in a 
S flavor state. 



Orbital angular momentum based on covariant 
derivatives 



Smearing of quark fields is based on the distillation 
method of Ref. pj)] . It is used in order to filter out the 
effects of small scale fluctuations of the gauge fields and 
it provides a spherically symmetric distribution of each 
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Dirac IR Emb 


p ® s 


S 1 1 

° 2 1 

1 1,2 


^75 ( + I/°)ms|s)ms + p)ma s)ma) 
|p)s|s)s 


M 1 1,2 

3 

3 1 

2 - 1 


MS |p)s|s)ms 

^5 ( — |p)ms|s)ms + |p)ma|s)ma) 
p)ms|s)s 


MA |p)s|s)ma 

^75 ( + |p)ms|s)ma + |p)ma|s)ms) 
p)ma «}s 


A 1 1 


( - p)ms s)ma + p)ma s)ms) 



TABLE VI: Symmetries of Dirac spin states based on direct products of p-spin and s-spin states for three quarks. States 
labeled as |p)z refer to a p-spin state with symmetry H from Eqs. ( A4| A6 1 . Similarly, states labeled as \s)t refer to a s-spin 
state with symmetry T from Eqs. ( A4|A6 1. Direct products of the p-spin and s-spin states yield sums of three-quark terms in 
which each quark has a p and s label. That determines each quark's Dirac index according to Table |V| 
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0MS 


0MA 
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1 

2 


+ k 





\udu) MS 


Md?i) MA 
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1 
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1 
2 





— dud) ms 


- \dud) MA 


As 








-1 


sud) M s - \uds) M s) 


swd) M A - uds) M A) 


^8 


1 


+ 1 


-1 


usu)ms 


|ms'u}ma 




1 





-1 


|wsd)MS 


MSd) MA 
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-1 


-1 


dsd) ms 


I dsd) MA 


■=■0 
—8 


1 

2 




-2 


— \SUS) MS 


— \SUS) MA 


"8 


1 
2 


i 

'?, 


-2 


- sds} MS 


| sds) MA 



Decuplet, 10 
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S 1 


4>s 


A++ 


3 
2 


+1 





\uuu)s 


A+ 


3 
2 


+ 2- 





\uud)s 
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3 
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1 
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\udd)s 


A" 


3 
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3 
2 





\ddd)s 


^10 
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+ 1 


-1 


\uus)s 


^10 
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-1 


|ttds)s 


^10 
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-1 


|dds)s 


=0 
"10 
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|sus)s 
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-2 


|sds)s 


0" 








-3 


|sss)s 



Singlet, 1 




I h 
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AS 





-1 


\dus) A 



TABLE VII: Flavor octet, decuplet and singlet constructions. 



quark field that carries no orbital angular momentum. 
In order to obtain higher spins, it is necessary to add 
spatial structure using covariant derivative operators, as 
described in the text, which are combined in definite sym- 
metries that correspond to orbital angular momenta in 
the continuum. For a single derivative, (D^), the two 
symmetry combinations are given in Eq. (J8|, while for 
two derivativesi, (D^), the combinations are given in 



Eq. (12) 



Operators that have good spin in the continuum are 
built by applying some number of derivatives to the Dirac 
spinors. Using the SU(2) Clebsch-Gordan coefficients 
to combine orbital and spin angular momenta, the one- 
derivative operators are, 

llM] = m ; S, s\ J, Af)4! i. m **'- 

171,8 

(A7) 

where S = | or | are the possible spin states of three 
quarks in the absence of derivatives. Reference [7] devel- 
oped single derivative operators for baryons and it pro- 
vides some examples of the incorporation of combinations 
of covariant derivatives into three-quark operators. 

Additional derivatives together with SU(2) Clebsch- 
Gordan coefficients are used to obtain higher J states. 
For example, the two-derivative operators are first com- 



bined to get L = 2 , 

D [ lU,m = E <!' TO i5 *> m ^ M >4Ml (A8) 



This L = 2 derivative operator is then applied to a spinor 
\E ,S ' S as follows, 



4i 2 * [s] 



J,M] 



m,s 



%m;S, S \J,M)D [ H 2 M * 



,S,s 



(A9) 

Derivative operator constructions for singlet, octet and 
decuplet follow from Table [VnTl The single-site operators 
are symmetric in space, flavor and Dirac indices. With 
one derivative, a mixed symmetry flavor and spin con- 
struction is combined with mixed symmetry derivative 
operators to make the overall symmetric combinations. 
Similarly, with two derivatives, various spin-flavor sym- 
metry states are combined with the derivative operators 
to make overall symmetric operators. 

As noted above, three quarks with no derivatives can 
form at most spin | and | states. The corresponding 
lattice irreps G\ and H are faithful representations of 
these continuum spins, and hence form the basis of the 
constructions of higher spins. The numbers of operators 
with up to two derivatives are shown in Table |IX| The 
number of operators classified according to total spin J 
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MA 0SXMA 


(10,2) 


l 


2 




4>sXms 






0SXMA 
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0MAXS 
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<"AmaXa 






^msXa 


(1,4) 




1 




</>aXma 






0aXms 


(1,2) 
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<^aXma 






0AXMS 


A (10,2) 
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</>SXA 




(8,4) 




1 




75 


</>msXma 


— 


/Wxms) 


(8,2) 
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75 


^msXma 
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TABLE VIII: Local operators classified according to symmetry of flavor and Dirac spin. The dimensionality of the SU(3) 
representation is shown. The number of Dirac spin embeddings (number of operators) in a non-relativistic (p = +) construction 
are shown in column A^ nonrc i, and the number of constructions featuring some number of lower components is shown in column 
iVfei. The total number of embeddings is the sum of N ncmTC i + N Te \. The multiplicity of operators in the non-relativistic case is 
56 (S), 70 (MS), 70 (MA) and 20 (A), and corresponds to the conventional non-relativistic 5*17(6) ® 0(3) construction. The 
relativistic construction, which involves both positive and negative parity operators, corresponds to the reduction of 5*17(12). 
Note that the flavor singlet operators are distinct from the octet and decuplet operators. In the 5(7(3) flavor limit, the flavor 
singlet states do not mix with the octet or decuplet states. 



and irrep are shown in Table |A 3| A general feature of 
the operator construction is that there is always an equal 
number of positive and negative parity operators. For ex- 
ample, for every operator in G\ g , there is a corresponding 
operator in G\ u and similarly for H g and H u . There are 
no single-site operators in G*2. 

These constructions provide operators that have good 
total angular momentum up to J = % in the continuum 
limit. However, they are reducible with respect to the 
octahedral group that represents the symmetry of a cubic 
lattice. 



elementary G\ and H octahedral irrep states, i.e., 



Appendix B: Quantum mechanics of continuum spin 
in the octahedral representation 



In this appendix we develop the subduction of SU (2) 
quantum states to irreducible octahedral states in the 
continuum. The lowest spins are trivial as suitable sub- 



ductions for spins J 



\ and J 



are provided by the 



Gx,t, 
H,r, 



r = to, 

2 



r = m. 

2 



(Bl) 



A suitable subduction for spin 1 also is trivial in terms 
of octahedral irrep T\, 

\[l,m]) = \T 1 ,r,[l}}, r = 2-m. 

Here and in the following we label octahedral states that 
carry continuum quantum numbers J, M by placing the 
quantum numbers in brackets, i.e., |[J,M]). The angu- 
lar momentum basis states are orthogonal in the contin- 
uum and consequently the octahedral states labelled as 
| [J, M]) are orthogonal in the continuum limit as follows, 



[J,M\ [J',M'])=5j t j,5 



JJ,J'OM,M'- 



(B2) 



However, they are reducible. 

Irreducible octahedral states, |A, r, [J]), that are sub- 
duced from a single spin J are labelled by the octahe- 
dral irrep, A, row r, and spin [J] in brackets. Examples 
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Singlet, 1 


IR 


d = 


d= 1 d = 2 


Total 


d 


1 


4 9 


14 


H 





5 17 


22 


G 2 





1 8 


9 


Octet, 8 


IR 


d = 


d = 1 d = 2 


Total 


Gi 


3 


8 17 


28 


H 


1 


11 36 


48 


G 2 





3 17 


20 


Decuplet, 10 


IR 


d = 


d= 1 d = 2 


Total 


Gi 


1 


4 10 


15 




2 


5 19 


26 


G 2 





1 10 


11 



TABLE IX: Numbers of singlet, octet and decuplet opera- 
tors for each parity according to the irrep and the number 
of derivatives, d. Total derivative constructions have been 
removed. 



Singlet, 1 


Rep 


J — ? J — 2 J — ^ J— 2 


Total 


Gi 
H 
G 2 


13 1 

13 8 1 
8 1 


14 
22 
9 


Octet, 8 


Rep 


J — 2 J — 2 J — 2 ^ ~ 2 


Total 


Gi 
H 
G 2 


24 4 

28 16 4 
16 4 


28 
48 
20 


Decuplet, 10 


Rep 


J— 2 J — 2 ^ ~ 2 ^ ~ 2 


Total 


Gx 
H 
G 2 


12 3 

15 8 3 
8 3 


15 
26 
11 



TABLE X: Number of singlet, octet and decuplet operators 
according to continuum spin and subduced irrep. With two 
derivatives, at most J = | can be reached. 



appear on the right side of Eq. (Bl). Owing to the or- 
thogonality of different octahedral irreps and rows, these 
states form an orthonormal set obeying, 



(A,r,[J] A',r',[J]) = 5 A<A ,5 r y 



(B3) 



A general octahedral irrep can contain an infinite number 
of continuum spins. States transforming as the same oc- 
tahedral irrep and row but subduced from different spins, 
such as \H,r, [§]), \H,r, [§]) and \H,r, [|]) 5 are distin- 
guished by their [J] labels. These states are orthogonal to 



one another as shown in Eq. (Bll ). The lowest spins have 
a one-to-one relation between the |[J, M)) and |A,r, [J]) 



labcllings as in Eq. (Bl ) but higher spins do not. 

Higher-spin states can be constructed from direct prod- 
ucts of lower-spin states by use of the SU(2) Clebsch- 
Gordan formula for direct products of states of spins Ji 
and J-i as follows, 

\[J,M})= \{Ji,m 1 })(x)\[J 27 m 2 }) 

771 1 ,7712 

x(Jimi; J 2 m 2 \JM), (B4) 

where ( J\m\\ J^m-i\JM) is a SU(2) Clebsch-Gordan co- 
efficient. Equation (B4) provides a block-diagonal uni- 
tary transformation between the basis of (2 Ji+1)(2 J2+I) 
product states on the right side and the equal number of 
states in the basis of total angular momentum on the left 
side for J in the range I Ji — J 2 | < J < J\ + J 2- 

Each of the octahedral states should be expanded in 
terms of a set of (2Ji + 1)(2J2 + 1) linearly independent 
states transforming as irreducible representations of the 
octahedral group and subduced from a single J value. 
Tables |Xl| and |XlT| show the relevant states based on using 
Ji = 1 and J 2 = I or | to construct J = | or J = | . 

The limited dimensions of the octahedral irreps require 
multiple occurrences of some irreps in the subduction of 
high spins. In the continuum, the different occurrences 
of the same irrep and row provide linearly independent 
states, as we will show by construction. Representa- 
tions with multiple occurrences are denoted as | ra A, r, [J]), 
where "A denotes the n th occurrence of irrep A, r denotes 
the row of the irrep and [J] shows the spin from which 
the state is subduced. When there is a single occurrence, 
the left superscript is omitted. 

Expanding in terms of a complete set of irreducible 
octahedral states gives, 



[J,M]) = J2 

"A,r 



n A,r,[J])( n A,r,[J}\ [J,M]), 



A, 



s 



J,M 
'A,r 



= ( n A,r,[J]\[J,M}), 



(B5) 



where we sum over all irreps, including multiple occur- 
rences of the same irrep, that are linearly independent. 
The subduction matrix, S^}^ r , is defined by the overlap 
of an irreducible octahedral state subduced from a single 
J and the spin state | [J, M] ) . Orthogonality proper- 
ties of the subduction matrices follow from substituting 
Eq. (lB5j) into Eq. ([Bl with J = J' as follows, 



EY, s ir^(^r,[j]Wy,[j]) = 5 



M.M' • 



A,r A' 



Because of Eq. (|B3j), this reduces to 

E C J,M C J,M' r 



M' 



(B6) 
(B7) 
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TABLE XI: Occurrences of octahedral irrep states with spins J = i, | and § based on using Ji = 1 and J2 
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7 
2 


IR states 




\H,r,\l]), \G 2 ,r,\l]) 


|Gi,r,m> 




. |G2,r,[|l> 



TABLE XII: Occurrences of octahedral irrep states with spins J = §, | and | based on Ji = 1 and J2 



Summation of the squares of S^ r over A and r gives the 
normalization condition, 



E C J,M C J,M 
^A.r ^A.r 



1. 



(B8) 



A.?- 



Because the | [J, M] ) states are a complete set over the 
subspace of spin J, there is a sum rule 



53<A,r, [J]|[J,M])([J,M]|A',r',[J]) 

= (A,r, [J]|A',/,[J]). 



M 



where Y^m \[JiM])(\J,M}\ = 1 was used to obtain the 
right side. The left side involves a sum over products 
of subduction matrices and the right side involves the 
8- functions of Eq. (|B3|) , yielding 



(B9) 



E C J,M C J,M c X 
6 A,r 6 A',r' = d A,A'<V 



Multiplying a subduction matrix times both sides of 
Eq. (|B5|) and summing over M yields 



M 

= I A, r 



a/ 



A',r' 



(BIO) 



where Eq. ( B9 ) was used to evaluate the summations on 
the right side. This last equation shows how the subduc- 
tion matrix is used. Once one has a realization of | [J, M]) 
states, the subduction matrix is applied to obtain the oc- 
tahedral irrep states that are subduced from a single J 
value. Using Eq. (BIO) leads to an important extension 
of Eq. pii 



(A, r, [J] |A',r', [J'] ) = J2 S J ^{ [J, M] I [J', M'] )sf,f, 



M,M' 

cJ,M x x cJ'M' 
s >->A .. Oj.j'Om.M' . , 



M,M' 



E C J,M C J,M r 
5 A,r ^AV' 5 



J, J' 



A/ 



= <5A,A"5r,r"5j : J' 



(Bll) 



where Eq. (|B9j) was used in the last step. The octahedral 
states subduced from single J values are orthonormal 



with respect to J as well as with respect to octahedral 
irrep, A, and row, r. 

Substituting the expansion in terms of orthonormal 
irrep states for each octahedral state in the Clebsch- 



Gordan formula of Eq. ( B4 ) gives 



El 

n A,r 



53 |A 1 ,ri,[J 1 ]>®|A 2 ,r 2 ,[J 2 ]) 



mi ,wi2 
Ai,ri 
A 2 ,i- 2 

x S J ^S J ^{j imx] J 2 m 2 \jM) 

(B12) 

Here the octahedral irrep states on the right side are 
assumed not to involve multiple occurrences so left su- 
perscripts are omitted for Ai and A 2 . 

The rules for combining direct products of octahedral 
group irreps are similar to those for continuum spins. The 
direct product T\ ® H corresponds to a direct product of 
J = 1 and J = I continuum irreps. It yields a G\ irrep 
that corresponds to spin-|, a H irrep that corresponds 
to spin- 1 and a pair of irreps, H' and G 2 , that taken 
together correspond to spin-[|]. The H and H' irrep 
states are orthogonal in the continuum limit because they 
are subduced from different spins, i.e., J = I and |. 

Direct products of octahedral states are equal to sums 
of irreducible states according to the Clebsch-Gordan for- 
mula for the octahedral group, 




|A l ,n,[Ji]>(8)|A 2 ,r 2 , [J 2 ]) = 

53 |AA!®A 2 ,r) • C 

Aa^Aj ,r 

(B13) 

where the states produced on the right side are labeled by 
the irreps involved in the direct products. The [Ji] and 
[J 2 ] labels of the states on the left side are passive and do 
not affect the direct products. This gives the following 
expansion, 



E 

"A,r 



"A, r, [J])S. 



J,M 
'A,r 



E 



|A Al ®A a ,r)ft 



J,M 



(B14) 
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TABLE XIII: Occurrences of octahedral irrep states in the subduction of spin | based on Ji = 1 and J% = | 



where matrix 72. 



J,M 

Aa 1 ® A 2 



r is defined by 



Using Eq. (|B5[), this becomes 



n 



J, At 



^ | Ai A 2 A Al!glA2 

Ai.ri 
A 2 ,r 2 

cJi,mi cJ 2 ,m 2 
X °A 1: n °A 2 ,r 2 

x ( Jimi; J2W2I JM). 



n A,r, [J]) = E - 4 AA 1 ®A 2 ,"A[j]|A Al( gA 2 ,r), 



where 



(B15) 



|AA 1( g.A 2 , r) 



The notation used here is based on the fact that the 
states produced by the Clebsch-Gordan expansion are 
general octahedral irrep states, not states subduce d from 
a single J value. For the construction based on Eq. (B12 1, 
the SU(2) spins in the range | Ji — J2I < J < Ji + Ji 
provide a complete set of {2J\ + 1)(2J 2 + 1) linearly in- 
dependent states and the octahedral irrep states on the 
right side of Eq. (B14) can be expanded as a linear com- 



A Al8 A 2 ,n 

At 



[J,M]). 



(B20) 



(B21) 



Table IXIVI shows the subduction matrices for the ele- 
mentary states of spin s , 1 and | . They are unit matrices 
in the basis used. Starting with the known subduction 
matrices for spins J\ = 1 and J 2 = §, the subduction 
matrix for spin | can be obtained as follows. Evalu- 



bination of the octahedral irrep states subduced from a 
single spin as follows, 



ate Eq. (B15) to obtain matrices TZ^ 



7? 5 

,r> H Tl ®H,r 



and 1Z^' M 



Matrix TZ 



H 



overall constant factor A 



.H 



is equal within an 

to subduction ma- 

M 



|A Al( g>A 2 ,r) 



— M — A'f - 1 1 

trix Si,\ . Similarly, 1Z^ T s>H<r and Kg7ti <*n ,r are ec l ual 



Hr . „„„„ "J! -~M Tl9H ,r 'J2T 1 

-At -At 

within overall factors to S^' r and Sq r 
all factors lead to S matrices normalized as in Eq. ( B8 1 



The over- 



J1+J2 

E "^A All8A2 , n A[j'] |'A, r, [J'] ) 
J'=\Ji-M 

(B16) 

A (g>A 2 an d "A [J'] The subduction matrices that are obtained for spin | are 

indices 

Substituting Eqs. ( |B13| > and ( |B16[ ) into Eq. ( |B12| ) and 
using Eq. (B5| gives 



where -4a Ai8A9 ,«a[,/'], is a matrix in the A 



given in Table |XV} The overall factors determine matrix 
"4"A[ j],a Ai , which is a unit matrix in this example, i.e., 
Eq. (B16l takes the form, 



J, At 
A,r 



"A,r 



E ^ 



" l A Al8A2 ,"A[J'] 

x|"A,r,[J']X« A2 , r (B17) 




1 N 
1 
,0 1, 



\H,r, [§]) 
l^.r, [|]> 
\|G 2 ,r, [f]>, 



Because the octahedral states | "A, r, [J] ) that are sub- 
duced from a single spin form an orthonormal set, their 
coefficients must be the same on both sides of Eq. (B17|. 



(B22) 



Thus, the subduction matrix for total spin, J, obeys 

-.j.m \ , r>JM 
2^ • A A Al8A2 ,'W[J]/c A 



A,®A 2 > r 



(B18) 



As will become evident, matrix A is orthogonal: ,A*A T = 
1. When the above equation is multiplied by \[J, M]) on 
both sides, and then summed over M, we get 

E«Pm]) = e 

At A Al g,j\ 2 



E<^ A2 ^|[ J ' M ])(B19) 

At 



This block-matrix equation holds when the row indices, 
r, are the same on both sides. Note that orthonormality 
of the states on the right side implies that the states on 
the left side are orthonormal. 

For higher spins, the Clebsch-Gordan coefficients of 
the octahedral group do not provide a block-diagonal 
result for the subduction, i.e., matrix A takes a non- 
trivial form. This is demonstrated for the construction 
of spin I based on J\ — 1 and J 2 = |, which involves 
use of the previous subduction, [|] — > H © G 2 . Three 
orthogonal occurrences of irrep H are produced in the 
Ti® [§] = Ti<S>(H®G 2 ) direct products and two orthogo- 
nal occurrences of G 2 are produced. Taking into account 
the fact that two different H irreps are produced by each 
T X ®H product, i.e., T X ®H -> Gi®H®H'®G 2 , the three 
H states that occur are Ht^h, H't^h ana ^ -^7i®g 2 - 



27 



The two different G 2 states that occur are G 2 t 1 <s,h and 

The procedure is similar to the spin | case. 
First one evaluates Eq. (B15) to obtain matrices 

,M 



K 



TV; 



,M 

! 2 T-l®H ,r 



, S ,M 



and TZ%' . Each TZ matrix corresponds to a single 

— M — Al — Al - 

S matrix of the set S^ i r Sfj' r , Sfj' r , Sfr r 



shows the linearly independent states for this case. A 
total of four orthogonal occurrences of H, two of G\ and 
two of G 2 arise from T x ® [|] = T t ® (d © H © G 2 ). 
Two of the H irreps reproduce the previous results for 
subduction of spins | and | . The other two H irreps are 
yet-to-be-determined subductions of spin | . The two G\ 
irreps provide subductions of spins | and |. The two 



• M <}*' M and 



G 2 irreps provide subductions of spins 



I and f. 



The 



new element for spin | is that two occurrences of irrep 



$G 2 ,r ■ In t hi s case the subduction matrices for spin £ can R &m unknown They are dete rmined as follows. First 



be obtained simply by imposing the normalization con 



dition 2^M S A ~ S * 



1. That gives the subduction 



matrices of Table |XVH 

The octahedral states that result from the Clebsch- 
Gordan formula are linear combinations of the states 
subduced from a single J value. This is expressed by 
the block-matrix equation, 



/|Gi Tl ®H,r}\ 
\H Tl ®H,r) 



\Ht 1 ^g 2 , t ) 
\G 2Tl ®H,r) 



where matrix ^4 Tl ®2 is 



/I 












A 1 



■2 X 



/\Gur, 

Iff^f]) 
Iff *•[§]) 
Iff *•[§]) 

<**r, [I]) 
V|G 2! r[|]>7 



(B23) 



n \ 
o 






(B24) 

Note that row indices, r, are the same on both sides of 



Eq. (B23 1 so that it connects the 18 states on the left side 



(considering all allowed values of the row indices) to 18 
states on the right side. The states on the right side are 
orthonormal as in Eq. (Bill and the rows of the matrix in 



Eq. ( B23 1 are orthogonal to one another. It follows that 



the states on the left side also form an orthonormal set 



by construction. Matrix equation (B23 1 takes the form of 



Eq. (B16) and it provides a nontrivial example of matrix 



A, demonstrating the relation between irreducible states 
resulting from direct products of octahedral states and 
the irreducible states that are subduced from a single 
spin. It is straightforward to solve for the irreducible 
states subduced from a single spin by applying A^ 1 = A T 
to both sides of Eq. (B23|. That demonstrates that the 



irreducible octahedral states subduced from a single J are 
linear combinations of the irreducible octahedral states 
resulting from the Clebsch-Gordan formula. 

The same reasoning can be applied to the subdu ction 
of spin | based on J\ — 1 and J 2 — \. 



Table XIII 



matrices TZ 



J,M 

Aa 1 ®a 2 > r 



are calculated. They are related to 

the subduction matrices by a matrix _4 Tl ®f , which also 
relates the states formed from the TZ and S matrices as 



as follows, 



in Eqs. (BIO) and (B21 ). That relation can be expressed 



= A 7 



\Gl T 1 «.Gi,r) N \ 

\G\ T!®H,r) 

\H\ TiigiGi)'"! 

\H Tl0H ,r) 

|^Ti«iH' r ) 
l#T 1 «>G 2 ,7') 

\G 2 t 1 ®h , r) 
\\G 2Tl ®G 2 ,r) / 



where matrix A Tl ®2 is 



'2 X 



(\Gur, [§]>\ 

Gi,r,[i]> 
'H,r, [§]) 
2 H,r, [§]) 
\H,r[i\) 
Iff *"[§]) 
G 2 ,r,[r\) 
V|G 2 ,r[|])7 



/ /20 




(B25) 



f J 



(B26) 



Constants a n and b n express the unknown parts of 
the matrix that connect the four H irreps on the left 
side to linear combinations involving | H, r, [|] ) and 

fff^f])- 

The upper-left block of the matrix equation can be 
solved for 



Gi,r, [f] 



G\T x ®Hif\ 



(B27) 



,A1 



Subduction matrix Sq r is determined by a linear com- 
bination of TZ matrices with the same coefficients as in 



28 



Eq. (B27), i.e., 



20 2^ M 
27 /v 'Git 1 ®g 1 ,'- 



7 > M 

—72. 

27 Git^ 



(B28) 



The middle block gives four equations for the H irreps. 
They can be reduced to three by making linear combi- 
nations of pairs of equations to eliminate the \H,r, [|]) 
terms, and then further reduced to two equations by mak- 
ing linear combinations of pairs of equations to eliminate 
the \H,r, [|]) terms. The resulting equations provide 
candidates for the 1 H and 2 H irrep states subduced from 
spin | as follows, 




Eq.(B30l give the values of 61,62,63 and 64. They com- 
plete the determination of matrix A and can be used to 
show that it is orthogonal. Subduction matrices are de- 
termined as linear combinations of 1Z matrices with the 
same coefficients that appear above for the 1 H and 2 H 
states, 



S 



2 



M 




100 ,X -H^ H ,r 



225 '^H TlSlH ,r 



49 7-) 2- M 
675 ,X -H TlfSH ,r 

T 1 ®G 2 ' T 



(B30) 



The resulting subduction matrices for spin | are given in 
TableEVlIl 



l H,r, [§] 



where a n and b n are combinations of the unknown con- 
stants a n and b n . There remain four unknown constants 
here and three equations that constrain them in order 
that states \ l H,r, [§]) and \ 2 H,r, [§]> are orthonormal. 
A one-parameter family of solutions exists. It is sufficient 
for our purpose to obtain a single solution by choosing 
6 X = 0. From the first equation we find 



TABLE XIV: Elementary subduction matrices Sq r , <S T ^. 



and S^' r 
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m 1 +10-1 
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1 



'H,r, [§] 



Hn 



iH,r 



H'rp 



(B29) 



where constant a 1 was determined by normalizing the 
state. The second equation is then used to obtain a state 
that is orthonormal to the first one, which determines 
constants a 2 and b 2 , and yields 



2 #^,[§])= yf£g\H TlfsGl ,r) + 

H T 1 ®H-> r ) + V § H Tl ®G 2 ,r 



The coefficients appearing in Eq. (B29) give the values 



of ai, 02,03 ( 04 = 0) and the coefficients appearing in 



TABLE XV: Subduction matrices S^' r and S^ r 
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TABLE XVI: Subduction matrices S£™ r , S£™ and S£ 2 
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TABLE XVII: Subduction matrices S£™, 5^™ and S 2 2 H 
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